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We present a detailed study of the combinatorial interpretation of matrix integrals, 
including the examples of tessellations of arbitrary genera, and loop models on random 
surfaces. After reviewing their methods of solution, we apply these to the study of var- 
ious folding problems arising from physics, including: the meander (or polymer folding) 
problem "enumeration of all topologically inequivalent closed non-intersecting plane curves 
intersecting a line through a given number of points" and a fluid membrane folding prob- 
lem reformulated as that of "enumerating all vertex-tricolored triangulations of arbitrary 
genus, with given numbers of vertices of either color" . 
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1. Introduction 

Our first aim in these notes is to convince the reader that matrix integrals, exactly 
calculable or not, can always be interpreted in some sort of combinatorial way as generating 
functions for decorated graphs of given genus, with possibly specified vertex and/or face 
valencies. We show this by expressing pictorially the processes involved in computing 
Gaussian integrals over matrices, what physicists call generically Feynman rules. These 
matrix diagrammatic techniques have been first developed in the context of quantum 
chromodynamics in the limit of large number of colors (the size of the matrix) [|l| , and 
more recently in the context of two-dimensional quantum gravity, namely the coupling 
of two-dimensional statistical models (matter theories) to the fiuctuations of the two- 
dimensional space into surfaces of arbitrary topologies (gravity) 0. These toy models 
for non-critical string theory are a nice testing ground for physical ideas, and have led to 
many confirmations of continuum field-theoretical results in quantum gravity. The purely 
combinatorial aspect of these models has often been treated as side-result, and we believe 
it deserves more attention, especially in view of some spectacular results. Indeed, the non- 
critical string machinery allows one to relate critical properties (such as singularities of 
thermodynamic quantities) of the fiat space statistical models to those of the same models 
defined on random surfaces 0] . This has led for instance to recent progress in the study of 
random walks, by using the inverse relation to deduce fiat space results from gravitational 
ones 0. 

Once the connection is made between a combinatorial graph-related problem and some 
matrix model, we still have to compute the integral. Various powerful techniques have been 
developed, mainly in the context of various branches of physics, to compute those integrals: 
the original one is orthogonal polynomials [0, but only applies to "simple" models. More 
general is the saddle-point technique p|, that however only allows for computing these 
integrals in the limit of large size of the matrices. 

In these notes, we wish to present applications of the combinatorics of matrix models 
to some specific questions arising in physics having to do with folding. Folding problems 
arise in biology and physics in particular when considering polymers or membranes. An 
ideal polymer is a chain of say n identical constituents represented by segments (chemical 
bonds) attached to one another by their ends (atoms), around which they can rotate freely. 
Membranes are two-dimensional generalizations of polymers, i.e. reticulated networks 
made of vertices (atoms) linked by edges (chemical bonds), either in the form of a regular 
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lattice (tethered membranes) or in the form of networks with arbitrary vertex valencies 
(fluid membranes). Ideally, imposing that all bonds be rigid, the only possibility for a 
polymer or membrane to change its spatial conflguration is through folding, in which 
atoms (for polymers) or bonds (for membranes) serve as hinges. Quantities of interest 
for physicists are thermodynamic ones, characteristic of the systems when their size is 
large. Mathematically, these correspond to asymptotics of say the numbers of distinct 
configurations of folding of polymers or membranes of given length or area, when the 
latter tend to infinity. Both in the case of polymer and fiuid membrane folding, we will 
present matrix models allowing for the calculation of such thermodynamic properties. 

The polymer folding problem is better known in mathematics as the "meander prob- 
lem", namely that of enumerating all the topologically inequivalent configurations of a 
closed nonselfintersecting plane curve intersecting a line through a given number of points. 
This problem apparently first emerged in some work by Poincare in the beginning of the 
century, and reemerged in various areas of mathematics [§[ , from recreational math- 
ematics to the 16th Hilbert problem to computer science to the theory of knots and links. 



In physics, the formulation as a folding problem and the relation to matrix integrals [10 



have brought very different developments both algebraic [|T2[ [|T3[ [jT5| and numerical 
11] |1^. As a highly non-trivial outcome of our study, we will present the exact 



meander asymptotics recently derived in |jT^ 



The notes are organized as follows. In Part A, we give a detailed presentation of 
the combinatorial interpretation of Hermitian matrix integrals, as generating tools for 
fatgraphs (Sect. 2), and proceed to review their various methods of solution (Sect. 3). 
Part B is devoted to the application of this fact to the problem of enumeration of all 
distinct compact folding configurations of a closed or open polymer, namely the number 
of distinct ways to fold a (self-avoiding) chain of identical constituents onto itself: this 
is also known as the meander or semi-meander problem. After defining the problem and 
reviewing a few known results (Sect. 4), we first present an algebraic formulation of 
the counting problem within the framework of the Temperley-Lieb algebra, omnipresent 
in the integrable statistical models, as well as the theory of knots and links (Sect. 5). 
This connection produces remarkable results, like the exact expression for the meander 
determinant, a meander- related quantity. In Sect. 6, we connect the meander and related 
problems to a multimatrix integral, that allows for the exact determination of meandric 
configuration exponents, in particular through the abovementioned connection ||^] beetween 
fiat and curved space models. Finally, we study in Part C a matrix model for generating 
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foldable triangulations, a sort of two-dimensional generalization of the meander problem, 
but without the self-avoidance constraint. These triangulations form a simple model for 
fluid membranes, encountered in physics and biology. After posing the problem (Sect. 7), 
we introduce a two-matrix integral generating these triangulations, that we solve in various 
ways (Sect. 8). Finally, Sect. 9 gathers a few concluding remarks. 



PART A: Matrix Model Combinatorics 



The aim of this introductory part is to familiarize the reader with the use of matrix 
integrals as tools for generating and enumerating graphs with various decorations. These 
correspond in turn to physical models of matter coupled to two-dimensional quantum 
gravity, in the form of fluctuating surfaces of arbitrary genus. In the following, we first 
present the combinatorial tools and give recipees to construct ad- hoc matrix integrals for 
various graph enumeration problems. We then expose various methods of computation of 
these matrix integrals, concentrating in particular on the planar graph limit. 



2. Matrix Models and Combinatorics 
2.1. Gaussian Integrals and Wick's Theorem 
Consider the following Gaussian average 



2.._ 1 r.-4 2n._.o. .^^_(2n)! 



Among the many ways to compute this integral, let us choose the so-called source integral 
method, namely define the Gaussian source integral 

1 r°° 2 2 
E(s) = (e^^) = ^ / e-^+'^da; = e^ (2.2) 



'2n 



oo 



Then the average ( |2.1| ) is obtained by taking 2n derivatives of S(s) = wrt s and 
by setting s = in the end. It is then immediate to see that these derivatives must be 
taken by pairs, in which one derivative acts on the exponential and the other one on the 
prefactor s. Parallelly, we note that {2n — 1)!! = {2n — 1)(2?7. — 3). ..3.1 is the total number 



Fig. 1: A star-diagram with one vertex and 2n out-coming half-edges stands 
for the integrand x^". In the second diagram, we have represented one non- 
zero contribution to (a;^") obtained by taking derivatives of by pairs 
represented as the corresponding connections of half-edges (the x's and s's are 
dual to one another, hence derivatives wrt s are in one to one correspondence 
with insertions of x in the integrand). 

of distinct associations of 2n objects into n pairs. We may therefore formulate pictorially 
the computation of ( ^.11 ) as follows. 

We first draw a star-graph (see Fig.l), with one central vertex and 2?! outcoming half- 
edges labelled 1 to 2n clockwise, one for each x in the integrand (this amounts to labelling 
the x's in a;^" from 1 to 2n). Now the pairs of derivatives taken on the source integral are 
in one-to-one correspondence with pairs of edges in the pictorial representation. Moreover, 
to get a non-zero contribution to {x'^"'), we must saturate the set of 2n legs by taking 
n pairs of them. Let us represent each such saturation by connecting the corresponding 
edges as in Fig.l. We get exactly {2n — 1)!! distinct closed star-graphs with one vertex. 
We may therefore write the one-dimensional version of Wick's theorem 



where the sum extends over all pairings saturating the 2n half-edges, and the weight 
is simply the product over all the edges formed of the corresponding averages (x^) = 
{(P / ds'^)'E{s)\s=o = 1. Each saturation forms a "Feynman diagram" of the Gaussian 
average. The edge pairings are called propagators (with value 1 here) . The Feynman rules 
are simply the set of values of these propagators. This may appear like a complicated 
way of writing a rather trivial result, but it suits our purposes for generalization to matrix 
models and graphs. Let us mention that the pictorial interpretation we have given for the 
computation of (x^"^) is not unique. For instance, we could have arbitrarily split the average 
into {x^x'^) for some integrers p, q with p + q = 2n. Then we would have rather represented 
two vertices with respectively p and q out-coming half-edges. Wick's theorem ( |2.3| ) states 
that we must now sum over (possibly disconnected) graphs obtained by saturating the 
p + q half-edges by pairs. The result of course is still the same. 




(2.3) 



pairings 
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As a last remark, and to make the contact with graphs, we may consider for instance 
the formal series expansion 



^(51,52,...) = (e^^^^'") (2.4) 



giX 



in powers of sr^, (72, whose coefficients are computed using (|2.1|) . Thanks to the previous 
pictorial interpretation, we may compute a typical term in the series expansion, say the 
coefficient (IK^*)^') of Ilfi'i^V^j^ by drawing Vi star-diagrams with i half-edges, i = 
1, 2, and saturating the 2E = half-edges in all possible ways by forming E pairs. 

Hence computing ( p.4|) is reinterpreted as a generating function of graphs with specified 
numbers of vertices of given valencies. 

2.2. One-Hermitian Matrix Model: Discrete Random Surfaces 

Let us now repeat the calculations of the previous section with the following Gaussian 
Hermitian matrix average of an arbitrary function / 

(/(M)) = -i— [ dMe-''^^'^f{M) (2.5) 



Zo{N) 

where the integral extends over Hermitian NxN matrices, with the standard Haar measure 
dM = Yli dMii Wi^j dRe{Mij)dIm{Mij) , and the normalization factor Zo{N) is fixed by 
requiring that (1) = 1 for / = 1. Typically, we may take for / a monomial of the 
form /(M) = Y[{i j)ei ^"ij' ^ ^ finite set of pairs of indices. Note the presence of the 
normalization factor (=the size of the matrices) in the exponential. Note also the slight 
abuse of notation as we still denote averages with the same bracket sign as in previous 
section: we may simply include the case of the previous section as the particular case of 
integration over 1x1 Hermitian matrices (i.e. real numbers) here. 

Like in the one-dimensional case of the previous section, for a given Hermitian NxN 
matrix S, let us introduce the source integral 

E(^) = (e^^(^*^))=e^ (2.6) 

easily obtained by completing the square - N{SM + MS) = [M - NSf - N'^S'^ 
and performing the change of variable M' — M — NS. We can use ( |2.6|) to compute any 
average of the form 
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Note the interchange of the indices due to the trace Tr{MS) = Yli^ijSji- As before, 
derivatives wrt elements of 5" must go by pairs, one of which acts on the exponential and 
the other one on the 5" element thus created. In particular, a fact also obvious from the 
parity of the Gaussian, (|2.7|) vanishes unless there are an even number of matrix elements 
of M in the average. In the simplest case of two matrix elements, we have 



Hence the pairs of derivatives must be taken with respect to Sij and 5"^^ for some pair z, j 
of indices to yield a non-zero result. This leads naturally to the Matrix Wick's theorem: 

( n ^^^■) = E n ^M^jM^i) (2-9) 

{i,j)el pairings P (ij),(A;Z)6P 

where the sum extends over all pairings saturating the (pairs of) indices of M by pairs. 

We see that in general, due to the restrictions ( |2.8D many terms in ( |2.9| ) will vanish. 
Let us now give a pictorial interpretation for the non- vanishing contributions in (p.9| ) . We 
represent a matrix element Mij as a half-edge (with a marked end) made of a double-line, 
each of which is oriented in an opposite direction. We decide that the line pointing from 
the mark carries the index z, while the other one, pointing to the mark, carries the index 
j. This reads 

M,j ^ \ ^ ^ 

(2.10) 

The two-element result (|2.8| ) becomes simply the construction of an edge (with both ends 
marked) out of two half-edges Mij and M^i , but is non-zero only if the indices i and j are 
conserved along the oriented lines. This gives pictorially 

i 

(2.11) 

Similarly, an expression of the form Tr(M"^) will be represented as a star-diagram with 
one vertex connected to n double half-edges in such a way as to respect the identification 



of the various running indices, namely 




(a) 





(b) 



(2.12) 



Fig. 2: An example of planar (petal) diagram (a) and a non- planar one (b). 
Both diagrams have n = 2p = 12 half-edges, connected with p = 6 edges. 
The diagram (a) has p + 1 = 7 faces bordered by oriented loops, whereas (b) 
only has 3 of them. The Euler characteristic reads 2 — 2h = F — E + 1 {V = 1 
in both cases), and gives the genus h — ior (a), and h — 2 ior (b). 



As a first application of this diagrammatic interpretation of the Wick theorem ( |2.9|) , 
let us compute the large asymptotics of (Tr(M'^)). To compute (Tr(M"^)), we must 
first draw a star-diagram as in ( |2.12| ), then apply (|2.9| ) to express the result as a sum over 
the saturations of the star with edges connecting its outcoming half-edges by pairs. To 
get a non- zero result, we must clearly have n even, say n = 2p. Note then that there are 
{2p — 1)!! such pairings, allowing us to recover the result of the previous section by setting 
A^ = 1, and simply replacing all oriented double lines by unoriented single ones. But if 
instead we take A^ to be large, we see that only a fraction of these {2p — 1)11 pairings will 
contribute at leading order. Indeed, assume first we restrict the set of pairings to "planar" 
ones (see Fig. 2 (a)), namely such that the saturated star diagrams have a petal structure 
in which edges only connect pairs of half-edges of the form (zj), i < j and (/c/), k < I with 
either j < k or k < i < j < I or i < k < I < j {in the labelling of half-edges, we have taken 
cyclic boundary conditions, namely the labels 1 and n + 1 are identified). In other words. 
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the petals are either juxtaposed or included into one-another. We may compute the genus 
of the petal diagrams by noting that they form a tessellation of the sphere (=plane plus 
point at infinity). This tessellation has V = 1 vertex (the star), E = p edges, and F faces, 
including the "external" face containing the point at infinity. The planarity of the diagram 
simply expresses that its genus h vanishes, namely 

2-2h = 2 = F - E + V = F +1-P =^ F = p+1 (2.13) 

Such diagrams receive a total contribution l/N^ from the propagators (weight 1/N per 
connecting edge), but we still have to sum over the remaining matrix indices ji, ^2, ■■■^jp+i 
running over the p + 1 oriented loops we have created, which form the boundaries of the 
F = p + 1 faces. This gives an extra weight N per loop (or face) of the diagram, hence 
a total contribution of N^'^^. So all the petal diagrams contribute the same total factor 
N to (Tr(M")). Now any non-petal (i.e. non-planar, see Fig. 2 (b)) diagram must have 
at least two less oriented loops, hence contributes at most 1/N to (Tr(M")). Indeed, its 
Euler characteristic is negative or zero, hence it has F < E — V = p — 1 and it contributes 
at most for N^~p < 1 /N. So, to leading order in N, only the genus zero (petal) diagrams 
contribute. We simply have to count them. This is a standard problem in combinatorics: 
one may for instance derive a recursion relation for the number Cp of petal diagrams with 
2p half-edges, by fixing the left end of an edge (say at position 1), and summing over the 
positions of its right end (at positions 2j, j — 1, 2, ...,p), and noting that the petal thus 
formed may contain Cj-i distinct petal diagrams and be next to Cp-j distinct ones. This 
gives the recursion relation 



Cp 



solved by the Catalan numbers 



J2cj-iCp-j co = l (2.14) 



_(2p)!_ 



Finally, we get the one-matrix planar Gaussian average 

hm 4(TV(M-)) = |^^P '{/"^^P (2.16) 
AT^oo N { otherwise ^ ^ 

This exercise shows us what we have gained by considering N x N matrices rather than 
numbers: we have now a way of discriminating between the various genera of the graphs 
contributing to Gaussian averages. This fact will be fully exploited in the next example. 
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We would now like to present an application of (2^) with important physical and 
mathematical consequences. Let us apply the matrix Wick theorem ( ^.91 ) to the following 
generating function /(M) = g^p{N ^-^^ giTr{M^)/i), to be understood as a formal power 
series of the gi, i — 1,2, 3, 4, ... 



{Ng 



ni ,n2 ,. • • r^O i> 1 



(2.17) 



E n 



(iVg, 



E 



ni,n2,...>0 'i>l labelled fatgraphs T 

with TiA i — valent vertices 



^-E{T) ^F{T) 



by direct application (|2.9| ). 




Fig. 3: A typical connected fat graph F, corresponding to the average 
(Tr(M)3Tr(M2)2Tr(M3)Tr(M4)2Tr(M6)Tr(MS)). The graph was obtained 
by saturating the ten star-diagrams corresponding to the ten trace terms, 
namely with ni = 3 univalent vertices, n2 = 2 bi- valent ones, — 1 tri-valent 
one, n4 = 2 four-valent ones, uq = 1 six-valent one and ns = 1 eight-valent 
one, hence a total oi V = IQ vertices. This graph corresponds to one par- 
ticular Wick pairing for which we have drawn the E = IQ connecting edges, 
giving rise to F = 2 oriented loops bordering the faces of F. 



In (|2.17|) , we have first represented pictorially the integrand n^(Tr(M*))'^" as a succession 
of rii z- valent star diagrams like that of (|2.12| ), i = 1,2, .... Then we have summed over 
all possible saturations of all the marked half-edges of all these stars, thus forming (non- 
necessarily connected) ribbon or fatgraphs F with some labelling of their half-edges (see 



Fig. 3 for an example of connected fatgraph). In (|2.17| ), we have denoted by E{T) the total 
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number of edges of F, connecting half-edges by pairs, i.e. the number of propagators needed 
(yielding a factor each, from ( |2.8|) ). The number -F(r) is the total number of faces of 
r. The faces of F are indeed well-defined because F is a fatgraph, i.e. with edges made of 
doubly oriented parallel lines carrying the corresponding matrix indices i = 1,2, ...N: the 
oriented loops we have created by the pairing process are interpreted as face boundaries, 
in one-to-one correspondence with faces of F. But the traces of the various powers of M 
still have to be taken, which means all the indices running from 1 to have to be summed 
over all these loops. This results in the factor per face of F in ( |2.17| ). Finally, the sum 
extends over all (possibly disconnected) fatgraphs F with labelled half-edges. Each such 
labelled graph corresponds to exactly one Wick pairing of ( |2.y| ). Summing over all the 
possible labellings of a given un-labelled fatgraph F results in some partial cancellation 
of the symmetry prefactors Yl^ l/(z"*7ii!), which actually leaves us with the inverse of the 
order of the symmetry group of the un-labelled fatgraph F, denoted by l/|Awt(F)|. This 
gives the final form 

j^v{r)-E{r)+F{r) 
S lAumi n«.-'^' (2.18) 

fatgrapha ' \ / \ 

r ~ 

where ni(F) denotes the total number of z-valent vertices of F and V{r) = ^^(F) is the 
total number of vertices of F. To restrict the sum in ( |2.18| ) to only connected graphs, we 
simply have to formally expand the logarithm of Zn, resulting in the final identity 

F^(^7i,^72,...) = LogZA.(^7i,^2,...)= E \Aut{T)\ (2-19) 



connected 
fatgraphs F 



where we have identified the Euler characteristic x(F) = F — E + V = 2 — 2/i(F), where 
h{r) is the genus of F (number of handles). Eqn. (|2.19|) gives a clear geometrical meaning 
to the Gaussian average of our choice of f{M): it amounts to computing the generating 
function for fatgraphs of given genus and given vertex valencies. Such a fatgraph F is in 
turn dual to a tessellation F* of a Riemann surface of same genus, by means of rii z-valent 
polygonal tiles, z = 1, 2, .... 



The result (|2.19| ) is therefore a statistical sum over discretized random surfaces (the 
tessellations), that can be interpreted in physical terms as the free energy of 2D quantum 
gravity. The name free energy stands generically for the logarithm of the partition function 
Zn. It simply identifies the Gaussian matrix integral with integrand f{M) as a discrete 
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sum over configurations of tessellated surfaces of arbitrary genera, weighted by some ex- 
ponential factor. More precisely, imagine only Qs = g while all other gi^s vanish. Then 
( p.l9|) becomes a sum over fat graphs with cubic vertices (also called cj)^ fat graphs), dual to 
triangulations T of Riemann surfaces of arbitrary genera. Assuming these triangles have 
all unit area, then n3(r) = A(T) is simply the total area of the triangulation T. Hence 
( p.l9|) becomes 



A{T) ^2-2h{T) 

connected triangulations T 



and the summand g^N"^ "^^ = e ^'^ is nothing but the exponential of the discrete version 
of Einstein's action for General Relativity in 2 dimensions, which reads, for any surface S 



SE(A,Af\S)=A [ (fxJW\ + — I d^xJ\^\R 
Js 47r Js 

= AA{S)+f^{2-2h{S)) 



(2.21) 



where g is the metric of S, R its scalar curvature, and the two multiplicative constants 
are respectively the cosmological constant A and the Newton constant A/". In ( |2.21| ), 
we have identified the two invariants of S: its area A{S) and its Euler characteristic 
x{S) = 2 — 2h{S) (using the Gauss-Bonnet formula). The contact with (|2.20|) is made by 
setting g = and = . If we now include all gi^s in (p.l9| ) we simply get a more 
elaborate discretized model, in which we can keep track of the valencies of vertices of F 
(or tiles of the dual F*). 

Going back to the purely mathematical interpretation of (|2.19D, we start to feel how 



simple matrix integrals can be used as tools for generating all sorts of graphs whose duals 
tessellate surfaces of arbitrary given topology. The size N of the matrix relates to the 
genus, whereas the details of the integrand relate to the structure of vertices. An important 
remark is also that the large N limit of ( p.l9| ) extracts the genus zero contribution, namely 
that coming from planar graphs. So as a by-product, it will be possible to extract from 
asymptotics of matrix integrals for large size N some results on planar graphs. 

2.3. Multi-Hermitian Matrix Case 

The results of previous section can be easily generalized to multiple Gaussian integrals 
over several Hermitian matrices. More precisely, let Mi, M2, ... Mp denote p Hermitian 
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matrices of same size N x N, and Qa,b, a,b = 1, 2, ...,p the elements of a positive definite 
form Q. We consider the multiple Gaussian integrals of the form 

where by a slight abuse of notation we still denote averages by the same bracket sign as 
before: the one-Hermitian matrix case of the previous section corresponds simply to p = 1 
and Qi,i = 1. The averages ( p.22|) are computed by extending the source integral method 



of previous section: for some Hermitian matrices Si, Sp of size N x N, we define and 
compute the multi-source integral 

and use multiple derivatives thereof to compute any expression of the form ( |2.22| ), before 
taking Sa 0. As before, derivatives wrt elements of the 5"s must go by pairs to yield 
a non-zero result. For instance, in the case of two matrix elements of M's we find the 
propagators 

{{Ma)^J{Mt)kl) = ^SuSjkiQ-^)a,b (2.24) 

In general we will apply the multimatrix Wick theorem 

( n (Mahj)= E n {iMahiM,)ki) (2.25) 

(a,i,i)eJ pairings pairs 

^ ' p (aij),{bkl)eP 

expressing the multimatrix Gaussian average of any product of matrix elements of the M's 
as a sum over all pairings saturating the matrix half-edges, weighted by the corresponding 
value of the propagator ( |2.24| ). Note that half-edges must still be connected according 
to the rule ( |2.11| ), but that in addition, depending on the form of Q, some matrices may 
not be allowed to connect to one another (e.g. if {Q~^)ab = for some a and b, then 
(MaMb) = 0, and there cannot be any edge connecting a matrix with index a to one with 
index b). 

This gives us much freedom in "cooking up" multimatrix models to evaluate generating 
functions of graphs with specific decorations such as colorings, spin models, etc... This is 
expected to describe the coupling of matter systems (e.g. a spin model usually defined on 
a regular lattice) to 2D quantum gravity (by letting the lattice fluctuate into tessellations 
of arbitrary genera). 
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An important example is the so-called gravitational 0{n) model |19], say on cubic 
fatgraphs. Its regular lattice version |j2^ is as spin model defined on the honeycomb 
(hexagonal) lattice, with only cubic (trivalent) vertices. A spin configuration of the model 
is simply a map from the set of vertices of the lattice to the n-dimensional unit sphere, 
hence the name inherited from the obvious symmetry of the target space. For n = 1, 
this is the Ising model, in which the spin variable may only take values ±1. For n > 2 
this gives an infinite continuum set of configurations, over which we have to perform a 
statistical sum, henceforth an integration over all the spin values. With a suitable choice 
of nearest neighbor interaction between adjacent spins, this statistical sum was shown to 
reduce to the following very simple loop model on the honeycomb lattice. The partition 
function Zo{n) of the regular lattice model is expressed as a sum over configurations of 
closed non-intersecting loops drawn on the honeycomb lattice, each weighed by the factor n, 
whereas each edge of each loop receives a Boltzmann factor K corresponding to the inverse 
temperature of the statistical mechanical model. The gravitational version thereof is very 
simple to figure out: we must replace the honeycomb lattice by its spatial fluctuations 
into arbitrary cubic fatgraphs. Then, on each fatgraph, we must consider arbitrary closed 
non-intersecting loop conflgurations, with the same weights as in the regular lattice case. 
The free energy of this model reads 

connected ' \ / I loop 

cubic fatgraphs F configurations I 

where as before we have inserted a weight g per cubic vertex (A(r) is the total area of 
the triangulation dual to F, also the total number of vertices of F), and for each loop 
configuration i on F, we have denoted by C{i) and E{i) respectively the total number of 
connected components of £ (number of loops) and the total number of edges occupied by 
the loops. Let us now present a multi-matrix model for ( p.26|) . Loop models can be easily 



represented by using the so-called replica trick: as we wish to attach a weight n per loop, 
we may introduce a color variable c = l,2,...,7i by which we may paint each connected 
component of loop independently. Summing over all possible colorings will yield the desired 
weight. Hence we are allowing the loops to be replicated, namely we may now choose among 
n different types of loops to build our configurations. To represent these configurations, we 
need n + 1 Hermitian matrices of same size N x N, say Ai,A2,..., An for the n colored loop 
half-edges, and B for the unoccupied half-edges. As loops are non-intersecting the only 
allowed vertices are of the form (unoccupied vertex) and A^B (loop of color c going 
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through a vertex, with exactly one unoccupied edge). The propagators must force loops 
to conserve the same color on each connected component, while unoccupied half-edges can 
only be connected to unoccupied ones, namely 

K 

{{Ac)ij{Ad)ki) = j^ScdSiiSjk 

{BijBki) = —Siidjk (2.27) 
{{A,),jBki) = 

where we also have added the weight K per edge of loop. This defines the inverse of the 
quadratic form Q we need for the corresponding matrix model, according to ( p. 24 ). This 
leads to the multi-matrix integral 



/rfSdAi...rfA,e-^Tr(V(B,Ai,...,A„)) 



fdBdAi...dAr,e-^T'<Vo{B,A„...,A^)) 

V{B, Ai, A.) = \B^ + i^Y.^l- - 9BC£ Al) 



2 2K 

c=l c=l 



(2.28) 



1 1 

Vo(S,Ai,...,AJ = -S2 + — 



2 2K 

c=l 



To get the free energy of the gravitational 0(n) model ( |2.26|) , we simply have to take 
the logarithm of ( |2.28| ), to extract the contribution from connected fat graphs only. As 
before, we may decorate the model with higher order terms of the form QiB'^/i to allow 
for arbitrary unoccupied vertex configurations and keep track of the fine structure of the 
fatgraphs, without altering the loop configurations. Another possibility is to remove the 
B^ term from ( |2.28| ), to restrict the loop configurations to compact ones, i.e. covering all 
the vertices of F. 

Another standard example, with a more complicated quadratic form, is the so-called 
gravitational g-states Potts model [lOfj defined say on arbitrary cubic diagrams F. Imagine 
that the diagram F is decorated by maps a from its set of vertices f (F) to 2^. We may 
alternatively view this as "spin" (or color) variables a G living on the vertices of F. The 
model is then defined by attaching Boltzmann factors corresponding to different energies 
of configurations of neighboring spins according to whether they are identical or distinct. 
More precisely, if v, v' are two vertices connected by an edge e in F, we define the weight 



w^{a) = e-^'^-(-).-(-') (2.29) 
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for some positive real parameter K (inverse temperature). The gravitational Potts model 
free energy is then defined as 

A{T) /\r2-2/i(r) 

F,ig\K,q)= Y: ' \Aut(T)\ ^ n-e(-) (2-30) 

connected ' \ / I maps edges 

cubic fatgraphs T CT:u(r)— >Zq e of T 

Let us now "cook-up" a multi- matrix model integral for ( p. 301) . We wish to generate 
all connected cubic fatgraphs F: this would be easily done by considering the example 
of the previous section with all g^s equal to zero except = g. To represent the spin 
configurations however, we need to be able to distinguish between q different types of 
vertices, according to the value of cr(f) G 2^. This is done by introducing q Hermitian 
matrices Mi, Mq with the same size N x N, with a Gaussian potential of the form 
NTv{Y,MaQabMb)/2, such that 

{Q-')ab = 1 + (e^ - l)5ab (2.31) 

in order for the propagators ( |2.24| ) to receive the extra weight when the "spins" a and 
b are identical, and 1 when they are distinct. Eqn.( ^^3l| ) is easily inverted into 

The free energy of the gravitational Potts model on cubic diagrams is therefore given by 

FNig\K,q) = LogZNig\K, q) 



ZN{g\K,q) j^Mi...dMqe-^^-Vo{M„...,M,) 
y(Mi,...,M,) = - J2 QabMM-^Y.^- 

a,b=l a=l 



(2.33) 



1 ^ 

Fo(Mi,...,M,) = 2 5Z Q-bMaMi 



a,b=l 



with Q as in ( |2.32| ). When q = 2 the model is nothing but the Ising model (also equivalent 
to the 0{n = 1) model, as mentioned above). To recover ( |2.3C1| ) from (|2.33| ), we simply 



apply the multi-matrix Wick theorem, and note that the overall contribution for each 
graph is ]Sf^~^+^ = jY2-2/i before, as we get a weight per vertex, 1/A^ per edge from 
propagators, and per oriented loop of indices, irrespectively of the spin values. 
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2.4- A generating Function for Fatgraphs 

In this section, we use the philosophy developed in the two previous ones to present 
a model of dually weighted fatgraphs, namely in which we wish to specify not only the 
structure of vertices but that of faces as well This uses the slightly more general 

notion of index-dependent propagators, obtained by considering one-Hermitian matrix 
models with Gaussian potentials of the form 

V{M) = \Y^ M.jQ^j.kiMki (2.34) 

ijkl 

where Q is some arbitrary tensor depending on the four matrix indices, but such that 
TrV{M) remains positive. Repeating the calculations of Sect. 2.3, we find propagators of 
the form 

{M,jM,i) = ^(^y^ + j^) (2.35) 

This gives even more freedom in tailoring specific models for combinatorial purposes. Let 
us now examine the choice 

Qij;ki = XiXj5ii5jk (2.36) 

for some Ai, Aa? > 0. We will also refer to the N x N diagonal matrix A with entries 
= Xidij. Then the propagator ( p.35| ) becomes 

{M.jMki) = -TT^SuSjk (2.37) 
iV AiAj 

If we consider the Gaussian average wrt ( |2.34D of the function f{M) = exp(—N Yli>i 9i^{M'^) /i) 
we find the following interpretation after taking the logarithm: 



i^iv(A;^i,^2,-) = Log(e ^«>i ^ ') 



N^-^^^^^ TT n.(r) TT V- 1 (2-38) 



^ \AutiT)\ J-J-^^ J- J- ^ y 

I ^ i>l faces i = l \ 



connected _ 
fatgraphs F / of F 



where, for each face / of F, we have denoted by m(/) the perimeter of its boundary 
(number of adjacent edges, or valency). The latter term comes from the summation over 
all the indices running over the oriented loops of the graph. But this time, we do not use 
^^j^ 1 = A^, but rather notice that each edge of a loop carrying the running index i comes 
with a factor 1/Xi from the propagators ( |2.37|) . Hence each face comes with a total factor 
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m{f) the face valency, which has to be summed over from z = 1 to A^, to yield 
( p.38|) . Introducing the variables 

N 



G, = Tv{A-'^) =Y,^ (2-39) 

^=l ^ 



we finally get the expansion 



F^iA;g,,g,,..)= Twf^ 11 ^r^'^^r^"^ (2-40) 

conn, fatgraphs F i>l 

where mi{r) denotes the number of z-valent faces of F. The expression ( |2.4C1| ) looks almost 
symmetric in gi and G^: it would actually be the case if the Gi were an infinite set of 
independent variables like the ^r^'s, but they are not. Due to the definition (|2.39| ), only 
the first G's are generically independent variables, all the others are dependent. Here is 
how we should read ( p.40|) . For any given A^, let us expand -Fat (A; (71, (72, •••) ^.s a sum over 
fatgraphs with < edges. Then the rhs only involves the numbers Gi, ...,Gn, as well 
as gi, ...,(7Ar, and is manifestly symmetric in the g^s and G"s, as the two sets of variables 
are exchanged by going to the dual graphs. The coefficient of each monomial corresponds 
to the graphs with specified genus, numbers of z-valent vertices, and numbers of z-valent 
faces. To get information on larger graphs, we simply have to take A^ larger and larger. 



3. Matrix Models: Solutions 

We hope that this series of examples has convinced the reader of the relevance of 

matrix integrals to graph combinatorics. Let us now see how to extract useful information 

from these matrix integrals. In this section, we will mainly cover the one-matrix integrals 

defined in Sect. 2.2. Multi-matrix techniques are very similar, and we will present them 

later, when we need them. More precisely, we will study the one- matrix integral 

. , /dMe-^Try(M) 

^^^^^ - /dMe-^TrVo(M) (3-1) 

2 . 2 

with an arbitrary polynomial potential V{x) = \ + Yli'iKiKk Vo{x) = This 

contains as a limiting case the partition function (|2.17[) of Sect. 2.2. Note also that we are 

not worrying at this point about convergence issueJi for these integrals, as they must be 

understood as formal tools allowing for computing well-defined coefficients in formal series 

expansions in the g^s. 

^ The reader may be more comfortable assuming Iie{gk) > to actually ensure convergence 
throughout this section. 
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3.1. Reduction to Eigenvalues 

The step zero in computing the integral ( |3.1j ) is the reduction to one-dimensional 
integrals, namely over the real eigenvalues mi, ...,mN of the Hermitian matrix M. This is 
done by performing the change of variables M — * (m, U), where m =diag(mi, m^r), and 
t/ is a unitary diagonalization matrix such that M = Umll^ , hence U G U{N)/U{1)^ as 
U may be multiplied by an arbitrary matrix of phases. The Jacobian of the transformation 
is readily found to be the squared Vandermonde determinant 

J = A{mf = Yl {mi- rujf (3.2) 

l<i<j<N 



Performing the change of variables in both the numerator and denominator of (|3.1|) we 
obtain 

J^N (imi...(imivA(m)^e~^^i=i ^^^^^ 



(3.3) 



3.2. Orthogonal Polynomials 

The standard technique of computation of ( |3.3| ) uses orthogonal polynomials. The idea 
is to disentangle the Vandermonde determinant squared interaction between the eigenval- 



ues. The solution is based on the following simple lemma: if Pm{x) = x"^ + XljLo Pmjx^ 
are monic polynomials of degree m, for m = 0, 1,...,A^— 1, then 

A(m) = det(mp^)i<ij<Ar = det(pj_i(mi))i<ij<Ar (3.4) 

easily derived by performing suitable linear combinations of columns. Let us now introduce 
the unique set of monic polynomials pm, of degree m = 0,l,...,A^— 1, that are orthogonal 
wrt the one-dimensional measure diJ,{x) = exp{—NV{x))dx, namely such that 



{Pm:Pn)= / Pm{x)Pn{x)dn{x) = hmSm,n (3.5) 

These allow us to rewrite the numerator of (|37^), using (|3.4| ), as 



N „ N-1 

e(ar)n / rf^(m,)p.(.)-i(mOp.(.)-i(m,)e-^^('"^) = iV! J] (3.6) 

o-.reSjv i=l i=0 
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We may apply the same recipee to compute the denominator, with the result A^! IljLo^ ^j^^ ' 
where the h^^"^ are the squared norms of the orthogonal polynomials wrt to the Gaussian 
measure dfxo{x) = exp{—Nx^/2)dx. Hence the /I's determine Zn{V) entirely through 



N-l , 

^^(^) =11 rk (3.7) 

i=0 

We may repeat this calculation in the presence of a "spectator" term of the form 
/(^i)/-^ some arbitrary function /. The term corresponding to z = 1 in this sum 
simply reads, by obvious modification of ( p^) 

<^r) [ d^(mi)^^^^p.(i)_i(mi)p.(i)_i(mi)e-^^("^^)x 
^ f 

X n / dn{m^)p^^i)_l{m^)pr(^)-lim,)e-^^^'^^^ (3.8) 

N-l N-l . .2 

^ n X ^ / .™/,„.) j: 

j=0 i=0 ' 

which is independent of the choice of index i {= 1 here). Therefore, the complete average 
wrt to the full matrix measure reads 



1^,, V d7m...dm^j^ Zl^ /(mOA(m)^e-^E ^(-0 

(3.9) 



N-l , .2 



Pi{my 



-NV{m) 



where we have added the subscript V to recall that the average is normalized wrt the full 
measure (using V instead of Vq in the denominator term). In the case where /(m) is a 
polynomial, this can be immediately rewritten in terms of the r's and therefore the /I's 
only, as we'll explain now. 

To further compute the /I's, let us introduce the two following operators Q and P, 
acting on the polynomials pm- 

QPmix) = XPm{x) 

o . ^ d , , (3-10) 

Ppm[x) = —Pm{x) 

dx 
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with the obvious commutation relation 



[p,g] = i (3.11) 

Using the self-adjointness of Q wrt the scalar product (/, (7) = / f{x)g(x)dfx{x), it is easy 
to prove that 

QPm{x) = XPm{x) = Pm+l{x) + SmPm{x) + rmPm-l{x) (3-12) 

for some constants and and that = if the potential V{x) is even. The same 
reasoning yields 

rm = , , m = l,2, ... (3.13) 

and we also set tq = ho for convenience. 

Moreover, expressing both {PprmPm) and {PprmPm-i) in two ways, using integration 
by parts, we easily get the master equations 

mi Pm-1) 



N {pm-liPm-l) (3.14) 

Q={V\Q)pmiPm) 



which amount to a recursive system for Sm and r^. Note that the second line of ( |3.14| ) is 



automatically satisfied if V is even: it vanishes as the integral over IR of an odd function. 
Indeed, this latter equation allows for computing the s's out of the r's, therefore becomes 
a tautology when all the s's vanish. Assuming for simplicity that V is even, the first 
equation of ( |3.14| ) gives a non-linear recursion relation for the r's. The degree k oi V 
actually determines the number of terms in the recursion, namely k — 1. So, we need to 
feed the k — 2 initial values of ro, ri, r2, Vks into the recursion relation, and we obtain 
the exact value of Zn{V) by substituting hi = rori...ri in both the numerator and the 
denominator of (|37?|). Note that for Vo(x) = x^/2 the recursion ( p.l4|) reduces simply to 

m _ {Qpl^\p^^l,) _ 

N (JO) JO) . 'rn 
\Pm-l-' Pm-l) 



and therefore h^m = h^Q^m\/N'^ = -\/27rm!/A^"^"'"^/^. The p^m are simply the Hermite 
polynomials. 
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Finally, the free energy of the model (|3.1| ) reads 



I N — l 

Fn{V) = LogZN{V) = NLogroJ ^ + V (iV - z)Log^ (3.16) 

i=l 



in terms of the r's. In the case when V is even and /(m) a polynomial, the average ( |3.9| ) 
is clearly expressible in terms of the r's only, by use of the recursion relation (|3.12| ) (with 
Sm = 0). Indeed, when /(m) = m", we get 



iV-l 



(lTr(M"))^ = 1 ^ ^^If^ (3.17) 
When n = 2 for instance, this simply reads 



N-l 



' =1 



3.3. Large N asymptotics I: Orthogonal Polynomials 

As mentioned before, the large limit of matrix integrals always has an interpretation 
as sum over planar (genus zero) fatgraphs. In particular, the free energy F]\[{V) of ( |3.16| ) 



can be expressed in an analogous way as ( p.l9| ) as a sum over connected fatgraphs, except 
that all but a finite number of g^s vanish, namely gm = for m > k + 1, and also for 
m < 2. The large N contribution is 

FN{V)^N^fo{V) + 0{l) (3.19) 

where foiV) is the planar free energy, namely that obtained by restricting oneself to genus 
zero fatgraphs. 

In view of the expression ( p.l6|) , it is straightforward to get asymptotics like ( |3.19|) , 
by first noting that as ho ~ \f^^ ^^^^ term in ( p.l6|) doesn't contribute to the leading 
order N"^ and then by approximating the sum by an integral of the form 

where we have assumed that the sequence tends to a function = r{i/N) of the variable 
z = i/N when N becomes large. This can actually be proved rigorously. The limiting 
function r{z) in (|3.20| ) is determined by the equations ( |3.14| ), that become polynomial in 
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this limit. In the case V even for instance, where V{x) = x"^ /'2 + J2^=2 92iX^^ / (2i) {k = 2p), 
we simply get 

i^) + i2C'7^)92Azy (3.21) 



z = r[ ^ 

i=2 



The function r{z) is the unique root of this polynomial equation that tends to z for small 
z (it can be expressed using the Lagrange inversion method for instance, as a formal power 
series of the (7's), and the free energy follows from ( |3.20| ). This allows also for computing 
the large limit of averages of the form ( ^.9| ) or ( |3.17[ ) . For instance, the large limit of 
( FT8|) reads 

1 

lim ( — Tr(M2))^ = 2 / dzr{z) (3.22) 

and more generally we have 







hm {-TriM-))v = { O lo d^^^)" ^ = 2p (3.33) 
N^oo^N ^ 1 otherwise ^ ' 



Note that in the case V = Vq (and r{z) = z) of Gaussian averages, (|3.23| ) reduces to the 
result (|2T|). 

3.4- large N asymptotics 11: Saddle- Point Approximation 

Let us now present another solution, which does not rely on the orthogonal polynomial 
technique nor requires its applicability. We start from the A^-dimensional integrals (|3.3|), 
that we rewrite 

^ ,,,, r rfmi...rfmAfe-^'^(""i'-'""^) 

where we have introduced the "actions" 

1 ^ 1 
S{mi,...,mN) = j^^V{mi) - — ^ Log 

1 ^ 1 

So{mi,...,mN) = —^Vo{mi) - — ^ Lo^ 



rrii — rrij \ 



(3.25) 



I rrii — ruj \ 



For large N the numerator and denominator of ( |3.24| ) are dominated by the semi-classical 
(or saddle-point) minimum of S and 5*0 respectively. For 5", the saddle-point equations 
read 

^ V'im,) = ^ Yl (3-26) 



dm-i N ^ — rrii 

J l<i<N J ' 
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for j = 1,2, ...,N. Introducing the discrete resolvent 

1 ^ 1 

--(-)-]v5:^ (3.27) 

i=l 

and multiplying (|]26D by l/(iV( z — rrij)) and summing over j, we easily get the equation 

1 ^ 



^ V'{mj)-V'{z) 



N ^ z — m-i 

j=i ^ 



- r 



AT2 



Assuming ljn tends to a differentiable function io{z) when — >^ oo we may neglect the 
last derivative term, and we are left with the quadratic equation 

uj{z)^ -V'{z)uj{z) +P{z) = 



P(z) 



1 ^ 
hm — > 



V'{z)-V'{m ,) (3.29) 

J = l 



N^oo N ^-^ z — rrij 



where P{z) is a polynomial of degree k — 2. The existence of the limiting resolvent uj{z) 
boils down to that of the limiting density of distribution of eigenvalues 

1 ^ 

p(^) = Jim^-^5(^-m,) (3.30) 

j=i 

normalized by the condition 



/ p{z)dz = 1 (3.31) 
Jr 



as there are exactly N eigenvalues on the real axis. This density is related to the resolvent 
through 



uj{z) 



f A±dx = ^ 4r / ^'^-'Pi^)dx (3.32) 
J z-x 2^"^ Jr 

m=l 

where the expansion holds in the large z limit, and the integral extends over the support 
of p, included in the real line. Conversely, the density is obtained from the resolvent by 
use of the discontinuity equation across its real support 

p{z) = lim Loiz + ie) — uiz — ie) z E supp(p) (3.33) 

2m e^o 
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Solving the quadratic equation ( |3.29|) as 

we must impose the large z behavior inherited from ( 3.31|)(|3.32|) , namely that oji^z^ ~ Xj z 
for large z. For k > 2, the polynomial in the square root has degree 2{k — 1): expanding 
the square root for large z up to order 1/z, all the terms cancel up to order with V'{z), 
and moreover the coefficient in front of 1/z must be 1 (this fixes the leading coefficient of 
P). The other coefficients of P are fixed by the higher moments of the measure p{x)dx. 
For instance, when k = 2 and F = Vb, we get 

Mz) = liz-Vz^-^) (3.35) 

It then follows from ( |3.33| ) that the density has the compact support [—2, 2] and has the 
celebrated "Wigner's semi-circle law" form 



Po{z) = —V'^-z^ (3.36) 

ZTT 

Viewing the resolvent ujq as the generating function for the moments of the measure whose 
density is po (through the expansion (|3.32|) ), we immediately get the values of the moments 

with Cp as in ( pj.l5| ). These are indeed immediately identified with the planar limit of the 
Gaussian Hermitian matrix averages (with potential Vo{x) = x'^/2) by using the following 
identity, valid for any V 

lim (^TrM'^)y = / x''p{x)dx (3.38) 

The result ( jO^ ) agrees with ( |2J^ ) and ([3:231) . Actually, comparing ( CT) with ( |3T7|) , 
we deduce the following expression for the limiting density p{x) in terms of the orthogonal 
polynomials of the previous solution 

p{x) = hm 1 y M^e-^^(^) (3.39) 

2 = * 

When V = Vq, we recover from ( |3.39p the Wigner's semi-circle from standard asymptotics 
of the Hermite polynomials. 
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In the general case, the density reads 



1 



p{z) = —^/AP{z)-{V'{z)Y (3.40) 

and may have a disconnected support, made of a union of intervals. It is however interesting 
to restrict oneself to the case when the support of p is made of a single real interval [a, 6]. 
It means that the polynomial V'{z)'^ — 4P{z) has single roots at z = a and z = b and that 
all other roots have even multiplicities. In other words, we may write 

iviz) = liV'iz) - Q{zW{z-a){z-h)) (3.41) 

where Q{z) is a polynomial of degree k — 2, entirely fixed in terms of V by the asymptotics 
uj{z) ~ 1/z for large 1^1 . For instance, for an even quartic potential, we have 

z^ z^ 

^(z) = -{z - gz^ - (1 - g— - gz'^)Vz^ - a^) 

2 (3-42) 
a^ = —(l-y/l-12g) 

p{z) = —(1 -g— - gz^)^o? - z2 

ZTT Z 



The planar free energy (|3.19|) is finally obtained by substituting the limiting densities 



p, po in the saddle point actions, namely 

/o(y) = ^(po,K))-^(p, V) 

f f , (3-43) 

S{p,V) = / dxp{x)V{x) — / dxdy p{x) p{y)Log\x — y\ 



This expression seems more involved than our previous result (|3.20|) , but is equivalent 
to it. In the case of the quartic potential of (|3.42|) , we find the genus zero free energy 
fo{g) ^ fo{V) 

foig) = ^Log ^ + ^(a^ - 4){a' - 36) (3.44) 



with as in (|]42]) 



In cases where the orthogonal polynomial technique does not apply (like in complicated 
multi-matrix integrals), however, the saddle-point technique always gives acces to the 
planar limit. 
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For completeness let us quickly mention the corresponding results for the gravitational 
0(n) model introduced in Sect. 2.3 (see |j2^ for details and a more general solution). We 
have to evaluate the large asymptotics of the partition function 



JdAi...dA^dBe-^^'(^o^B,A,,A2,-,A„)) 



W{B, A^, A2, A^) = V{B) + 1 ^ - ^ A 

^ c=l 

1 1 " 
Wo{B, Ai, A2, A^) = + ^ 5^ ^ 



(3.45) 



c=l c=l 

2 



c=l 



a simple generalization of ( 2.28 ) with some arbitrary potential V{B). Note that W is 



simply quadratic in the A^, so we can perform the Gaussian integrals overs the A's first. 
More precisely, the potential takes the form 

^ n N 

TrM^(S,Ai,A2,...,A^)=Try(S) + -5^ ^ {A,\kCl^k■A^ 



c)lj 

c=l ijkl = l 



Qik;ji = Sij5ki — g{5ijBki + SkiBij 
More compactly, the quadratic form reads 

Q = I (g) I ~ g{I (g) B + B I) (3.47) 

The Gaussian integral over the A's, normalized like in ( |3.45| ) simply gives det(Q)~^/^ for 
each integral, hence 

Z= ^ [ dBe-^^'^^^^det{l0l-g{I^B + B0l)y/^ (3.48) 

So the 0(n) model reduces to a one-matrix model, but with a complicated potential. We 
may now apply both the reduction to an eigenvalue integral and the large technique 
sketched above. We get the eigenvalue integral 

1 ^ 1 
S{h,,...M) = ^Y.^{h)-—^ J2 Log|6,-6,| 



+ ^ 5Z ^og{l-g{h + hj)) 

l<i,j<N 
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(3.49) 



Note that in the last term we may impose the constraint i ^ j, as the terms i = j only 
contribute for 0{jj) to S. The corresponding saddle-point equations db^S = read 

for i — 1,2, A^. Setting Xi — 1 — 2gbi, and v'{x) = j^V'{b), we get 

Multiplying this by j^i j^^ — Siud summing over i, we arrive at the large quadratic 
equation: 

-Q{z) + v'{z)uj{z) + v'{-z)uj{-z) = uj{z)'^ + uj{-z)'^ + nuj{z)uj{-z) 

1 ^ 1 

uj{z) = lim — 

N^oo N f-'^ z - Xi (3.52) 

^ v'(z)—v'(xj) v'(—z) — v'(Xj) 



Q{z) = lim 1 y 



N^oo N ^ Z — Xi Z -\- Xi 

1=1 

while the saddle-point equation turns into a discontinuity equation across the support of 
the limiting eigenvalue distribution p: 

uj{z + iQ) + uj{z — iQ) + nuj{—z) = v' [z) z G supp(p) (3.53) 

For generic n, u){z) is well defined only in an infinite covering of the complex plane, the 
transition from one sheet to another being governed by (|3.53|) . However, setting 



n = 2cos(7rz/) (3.54) 

we easily see that iiv = r / sis rational, this number becomes finite. Indeed, when expressed 
in terms of the shifted resolvent 



, , nv'i—z) — 2v'{z) , ^ , , 



the discontinuity and quadratic equations ( p.53|) and ( p.52| ) become respectively 



wi^z + zO) -|- w{z — zO) -f nw{—z) = 

(3.56) 

w{z)'^ -\- w{—zY + nw{z)w{—z) = P{z) 
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where 

P{z) = -^{nv'{z)v'{-z) - v'izf - v'{-zf) - Q{z) (3.57) 
— 4 

Introducing the two functions 

w±{z) = e^''''^w{z) + e^'^'iw{-z) (3.58) 

namely with W-{z) = w^{—z), we immediately get from ( p.56|) that 

w+{z + iO) = -e''^''w-{z - iO) 

w-{z + i0) = -e~''"'w+{z~iO) (3.59) 
w^{z)w-{z) — P{z) 

For rational u = r/s, we see that the combination 

w+{zy + {-ly+'w- {zY = 2S{z) (3.60) 

is regular across the support of p. For instance, if v{x) is polynomial (resp. meromorphic) , 
so are P{z) and S{z). We may now solve ( |3.59| ) ( ^^60|) for w± a-s 



w± 



(z) = iS{z) ± ^S{zy-P{zy)^ (3.61) 



and finally get the resolvent using (|3.58| ) ( |3.55|) : 

^^^^ _2v'{z) -nv'{-z) 1 ^ 

n2-4 2zsin(7rf) (362) 

X (e'^'^iSiz) + ^/S{zY - P{zy)^ - e-'^'^iSiz) - S{zY - P{zy)^) 

Assuming that v is polynomial, we must further fix the coefficients of P and S by imposing 
the asymptotic behavior u!{z) ~ 1/z for large \z\. If deg(f) = /c, then deg(P) = 2{k — 1) 
and deg(5') = s{k — 1), and P and 5" are entirely fixed if we impose moreover that p has a 
support made of a single interval [a, b] . 

3.5. Critical Behavior and Asymptotic Enumeration 

We have seen in Sect. 2 how to write the generating functions for various families of 
(decorated) fatgraphs in terms of Hermitian matrix integrals. It is a standard fact that 
the critical properties of these generating functions can be translated into asymptotics 
for these numbers of fatgraphs, say for large numbers of vertices. Indeed, writing such 
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a generating function as F{g) = J2n>o ^ri9^ , where denotes a number of connected 
(decorated) fatgraphs with n vertices, assume F has a critical singular part of the form 

F{g)s^ng--{9c-9f-^ (3.63) 

when g approaches some critical value gc, then for large n the numbers Fn behave as 

F„ ~ grn-'+'^ (3.64) 

Let us apply this to the example of (|3.1| ) with the quartic potential V{x) = \ — g\- 
The generating function Ff^i^g) = LogZ';v(^) = J2h>o f^io) decomposes into the 

generating functions Fh{g) 

fh{9) = J29''h,n (3.65) 

n>0 

where fh,n denotes the number of fatgraphs of genus h with n vertices. Using the genus zero 

answer for fo{g) ( |3.44| ), we find that the critical singularity is attained when g = gc = 

in which case f 0(9) sing ~ {9c — 9)^^"^ and 7 = —1/2 in ( |3.63D . This implies the following 

asymptotics for the numbers of connected genus zero fatgraphs 

12" 

/o,n T (3.66) 

n2 

More refined applications of the orthogonal polynomial techniques lead to the genus h 
result 

12" 

- ^i+lH-H) (3.67) 
and to a score of interesting behaviors for the numbers of decorated fatgraphs for specific 
potentials V, as well as for the 0(n) model with also specific potentials, for which any 
positive rational value of —7 in ( |3.63| ) may be reached. The exponent 7 is called the string 
susceptibility exponent. More generally, it has been shown that the coupling of a critical 
matter theory to gravity results in a behavior ( |3.63| ) of the genus zero free energy with 7 
given by the formula 

c-1- J(25-c)(l-c) 
7 = (3.68) 

where c is the central charge of the conformal field theory underlying the flat space critical 
matter model [^. The exponent 7 = —1/2 is characteristic of the "pure gravity" models, 
in which there is no matter theory, in other words c = 0. In the case of the so-called dense 
critical phase of the 0{n) model, one has 

e2 



and therefore 



c(n) = 1 — 6 , n = 2 cos(7re) (3.69) 

1 — e 



from (|3.68| ). This is confirmed by the saddle-point results ||2^ mentioned above. 
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3.6. Gaussian Words 



In the case of multi-matrix integrals, one may wonder how the beautifully simple 
result ( p.l6| ) ( ^^37|) is generalized. More precisely, we are interested in the large limit of 
the multi-Gaussian average of the trace of any word: 

rrfMi...rfMfce"^^'^^Sti ^^^')Tr(M7i...Ml^'=M"'^+\..M"'"M (3-71) 
= lim ^ — ^ 

N J dMi...dMke ^^i=i ^' 
In this setting, the result (|2.16| ) corresponds to k = m = 1 and reads 



[0 otherwise 

for all n > 0. The simplest way of computing ( |3.71| ) in a given situation is the blind 
application of Wick's theorem (|2.25|) , by further restricting the Wick pairings to the planar 
ones only, as — oo. For instance, for k = 2 and m = 2, we have 

Vn„U2,ns,n4 = ( (Tr (Mf ) ) ) 

(3.73) 

where we have isolated respectively the pairings in which elements of the blocks M^'^ and 
M^^ may be connected, then those in which elements of the blocks and M^"* may be 
connected, and we have subtracted the terms in which all pairings take place within the 
same powers of the M's to avoid counting them twice (once in each of the previous terms). 

To compute the most general trace of a word in say k Gaussian Hermitian matrices, 
we simply have to apply the following quadratic recursion relation, for lv = e^'^'^^^ 

mk — l 

Vni,...,nmk ~ ~ ^ ] ^ ^ni,...,ni^ni_|.i,...,n^fe (3.74) 

This relation is a compact rephrasing of the Wick theorem in the case of planar pairings, 
as the reader will check easily (the only relation to be used is X]o<i<fc-i ~ ^' E3 
for details). This gives a priori access to the average of any trace of word. The numbers rj 
are natural multi-dimensional generalizations of the Catalan numbers ( p.l5| ). 
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PART B: Polymer Folding and Meanders 

We now present a first important appfication of matrix integrals to the fundamental 
combinatorial problem of meander enumeration, also equivalent to the compact folding 
problem of polymers. After defining meanders and reviewing various approaches to their 
enumeration, we define a matrix model for meanders and discuss its solution in some 
particular cases. We finally present a general argument leading to the exact values of the 
meandric asymptotic configuration exponents. 

4. Folding Polymers: Meanders 

4.1. Definitions and Generalities 




compact . 
polymer meander 

folding 

Fig. 4: A typical polymer with 2n = 10 segments is depicted, together with 
one of its compact folding configurations, in which the segments have been 
deformed and slightly pulled apart for clarity. The corresponding meander is 
obtained by drawing an horizontal line (river) intersecting the folding configu- 
ration (road) through 2n = 10 points (bridges). The connected half segments 
of polymer have been replaced by semi-circular portions of road for simplicity. 

A polymer is a chain of say m identical constituents, modelled by rigid segments 
attached by their ends, that serve as hinges in the folding process. A useful representation 
to bear in mind is that of a strip of stamps of size 1 x m, attached by their ends. Given a 
closed polymer chain with m — 2n segments (i.e. in which the chain is closed so as to form 
a loop, see Fig. 4), we address the question of enumeration of all the ways of compactly 
folding the chain, in such a way that in the final (folded) state all the segements are piled 
up onto one another. 
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As illustrated in Fig. 4, we may keep track of the folding by intersecting the final 
state with an oriented line (say east to west), and slightly pulling apart the segments of 
the folded polymer. Redrawing the result with a horizontal east-west oriented line, and 
representing any two connected half-segments by semi-circles, we arrive at a meander, i.e. 
the configuration of a non-intersecting loop (road, made of semi-circular bits) crossing the 
horizontal line (river, fiowing east to west) through a given number m = 2n of points 
(bridges), as depicted in Fig. 4. The number of bridges will also be called the order of the 
meander. The number of distinct meanders of order 2n is denoted by !§]■ 

We extend this definition to a set of k roads namely meanders with k possibly in- 
terlocking connected components and a total number of bridges 2n (corresponding to the 
simultaneous folding of k possibly interlocking polymers with a total number of 2n seg- 
ments). The number of meanders of order 2n with k connected components is denoted by 

(k) 

. Note that necessarily 1 < k < n. These numbers are summarized in the meander 



= Y.^S (4.1) 



polynomial []TC 

n 

LVl 
k=l 

Given a meander of order 2n, the river cuts the plane into two parts. The upper part 
is called an arch configuration a of order 2n: it is a set of n non-intersecting semi-circles 
drawn in the upper half plane delimited by the river and joining the 2n equally spaced 
bridges by pairs. The lower part is clearly the refiection 6* of another arch configuration 
b of order 2n. Let denote the set of arch configurations of order 2n, and for any 
a,b E let c(a, 6) denote the number of connected components of road obtained by 

forming a (multi-component) meander with a as upper half and 6* as lower half. Then we 
have 

a,beA2n 

The number of distinct arch configurations is nothing but the Catalan number 

(2nV 

\A2n\=Cr.= . , (4.3) 

as A2n can be easily mapped onto the planar Wick pairings of a single star with 2n branches 
leading to (|2.16|) : indeed, we just have to pull all the bridges together so as to form a 2n- 



legged vertex; this results in a planar petal diagram like that of Fig. 2 (a). As an immediate 
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consequence we get m2n(l) 



4 /(yrn ) by use of Stirling's formula for 



large n. In general, the meander polynomial is expected to behave for large n as 



m2n{q) ~ c(g)- 



n 



a{q) 



(4.4) 



and we just proved that -R(l) = 4, a(l) = 3 and c(l) = I/tt. In Sect. 6.4 below, we 
will present an argument leading to the general determination of the exact value of the 
meander configuration exponent a{q) as a function of q. Note that 



s{q) = lim —Logm2n{q) = Log R{q) 



(4.5) 



is nothing but the thermodynamic entropy of folding per segment of a multi-component 
closed polymer. 






open 
polymer 



compact 
folding 



semi-meander 



Fig. 5: An open polymer with 8 segments, together with one typical com- 
pactly folded configuration, and its semi-meander version of order 9, obtained 
by intersecting the 8 segments of the configuration with a circle (represented 
in dashed line) , that also intersects the support (wall) to which the open poly- 
mer is attached. The free end of the folded polymer is then stretched into a 
straight half-line, that becomes the semi-meander's river, the free end becom- 
ing its source. The dashed circle is deformed into a winding road, crossing 
the river through 9 bridges. For simplicity, we have represented all the pieces 
of the road joining the various segments by means of semi-circles. Note that 
here the road winds three times around the source of the river. 



We may also address the problem of compactly folding an open polymer chain of n — 1 
segments [|10[ , attached to a wall by one of its ends (think of a strip of stamps, attached in 
a book by its left end). As illustrated in Fig. 5, we may intersect the folded polymer with 
a circle that also crosses the book's end. Extending the polymer itself into a half-line, and 
deforming the circle accordingly without creating new intersections, we arrive at a semi- 
meander of order n, i.e. the configuration of a closed non-self-intersecting loop (road, the 
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former circle) crossing a semi-infinite fine (river with a source, the former polymer+book) 
through n points (bridges). 

Note that, in a semi-meander, the road may wind around the source of the river, 
as illustrated in Fig. 5. We denote by M„ the number of topologically inequivalent semi- 
meanders of order n, and by M^*^^ the number of semi-meanders with k connected compo- 
nents, 1 < A; < n. We also have the semi-meander polynomial 

n 

fflM = Y^^nU^ (4.6) 



semi-meander meander 

Fig. 6: A semi-meander of order n is opened into a meander of order 2n. 
Think of the lower part of the river as pivoting around its source to form a 
straight line, while all bridge connections are deformed without any new in- 
tersections. The lower part of the meander is a rainbow, made of n concentric 
semi-circles. The winding number is just the number of semi-circular arches 
passing at the vertical of the center (former source). 

A semi-meander may be viewed as a particular kind of meander by opening the river 
as sketched in Fig. 6 so as to double the number of bridges n 2n, and by connecting 
them as they were in the semi-meander in the upper-half of the meander, and through 
the reflection of a rainbow arch conflguration in the lower one, made of n concentric 
semi-circles. The semi-meander polynomial is easily rewritten as 



fhr^iq) = ^47) 
aeA2n 



and we have the value at q — 1: m2n(l) = ~ A'^ /{^/^Tn^^'^). We also expect the 
asymptotics 



^^n{q) - c{q)^ (4.8) 



35 



and we have ^(1) = 4, a(l) = 3/2 and c{q) = l/^/n. Note again that the semi-meander 
configuration exponent a{q) wiU be determined as a function of q in Sect. 6.4 below. 

Conversely, a meander may be viewed as a semi-meander with no winding. It is 
therefore sufficient to solve the semi-meander enumeration problem, provided we keep 
track of the winding numbers. 



4-2. A Simple Algorithm; Numerical Results 



Multi-component semi-meanders of order n are in one-to-one correspondence with arch 



configurations of order 2n { \i.7i ). Our algorithm [|TT| is based on a recursive construction of 
arch configurations that allows to keep track of both the numbers of connected components 
and windings. To build an arch configuration of order 2n + 2 from one of order 2n, we 
may do either of the following transformations, both involving the addition of two bridges 
along the river, respectively to the left and right of the previous ones: 
(I) For each external arch (semi-circle contained in no other) connecting say the bridges 

i and j, replace it by two semi-circles, one connecting the new leftmost bridge to i, 

and one connecting j to the new rightmost bridge. 
(II) Add a large external arch connecting the two added bridges: it circles the whole 

previous arch configuration. 




Fig. 7: The tree of semi- meanders down to order n = 4. This tree is 
constructed by repeated applications of the transformations (I) and (II) on 
the semi- meander of order 1 (root). We have indicated by small vertical 
arrows the multiple choices for the process (I), each of which is indexed by 
its number. The number of connected components of a given semi-meander 
is equal to the number of processes (II) in the path going from the root to it, 
plus one (that of the root). 
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The corresponding semi-meanders are obtained by completing this arch with the re- 
flection of the rainbow r2n+2 as lower part. It is easy to show that applying (/) — (//) to 
all of yields exactly A2n+2- Moreover, (/) preserves the number of connected compo- 
nents, while (//) obviously increases it by 1 (the net result is to add a circle around the 
semi- meander). Applying successions of (/) — (//) on the "root" (semi-meander of order 
1), we may build the tree of semi- meanders of Fig. 7. 

The above algorithm is easily implemented on a computer and yields all semi-meander 
numbers with fixed winding and connected components up to some quite large orders. We 
give examples in Table I below. 



/ I 


M 


/ L 


M 


k 

rv 


1V121 


h 

rv 




1 


1 


16 


1053874 


1 


369192702554 


16 


2376167414 


2 


1 


17 


3328188 


2 


2266436498400 


17 


628492938 


3 


2 


18 


10274466 


3 


6454265995454 


18 


153966062 


4 


4 


19 


32786630 


4 


11409453277272 


19 


34735627 


5 


10 


20 


102511418 


5 


14161346139866 


20 


7159268 


6 


24 


21 


329903058 


6 


13266154255196 


21 


1333214 


7 


66 


22 


1042277722 


7 


9870806627980 


22 


220892 


8 


174 


23 


3377919260 


8 


6074897248976 


23 


31851 


9 


504 


24 


10765024432 


9 


3199508682588 


24 


3866 


10 


1406 


25 


35095839848 


10 


1483533803900 


25 


374 


11 


4210 


26 


112670468128 


11 


619231827340 


26 


26 


12 


12198 


27 


369192702554 


12 


236416286832 


27 


1 


13 


37378 


28 


1192724674590 


13 


83407238044 






14 


111278 


29 


3925446804750 


14 


27346198448 






15 


346846 






15 


8352021621 







Table I: The numbers ' of semi- meanders of order n with k connected components, 
obtained by exact enumeration on the computer: on the left, the one-component semi- 
meander numbers [k = 1) are given for n < 29; on the right, n is fixed to 27 and 1 < k < n. 
These numbers allow in turn for probing the asymptotics ( [4.4|)(|4^ , as functions of q. 
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In particular, we have the large q asymptotics of the radii 



R{q) = 2^{1 + ^ + ^ 



3 



29 
8g4 



81 



16g' 



R{q) 



2 

= g + l + - + 


2 2 
g^ g^ 


4 8 
g^ g6 


12 
"7 


10 4 12 


98 154 
+ qi2 + qi3 + 


124 10 

^14 + ^15 


102 

qi6 + 


20 




6 











+ 



46 
31 



(4.9) 




Fig. 8: The functions R{q) and -R(g) for < g < 4 as results of large n 
extrapolations. The two curves coincide for < g < 2 and split for q > 2 
with R{q) > R{q)- Apart from the exact value -R(l) = -R(l) =_4, we find the 
estimates R{0) = 3.50(1), R{2) = 4.44(1), R{3) = 4.93(1) and i?(4) = 5.65(1). 



At finite values of q, the numerical results displayed in Fig.8 reveal an interesting phase 
transition between a phase for q < qc of irrelevant winding, i.e. the numbers of meanders 

and semi-meanders are then asymptotically equivalent R{q) = -R(g)), and a strong winding 
phase q > q^ where the winding number is proportional to n, and therefore R{q) » R{q)- 
The transition point is estimated as gc — 2 with poor precision. We will propose an exact 
value for qc in Sect. 6.4 below. 

5. Algebraic Formulation: Temperley-Lieb Algebra 
5.1. Definition 

The Temperley-Lieb algebra of order n and parameter g, denoted by TLn{q), is defined 
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through its n generators 1, ei, 62, e„_i subject to the relations 

(i) ef = qci i = l,2,...,n- 1 

(ii) [ei,ej] = if |i- j| > 1 (5.1) 
{Hi) eiei±iei = z = 1, 2, n - 1 

This definition becomes clear in the "domino" pictorial representation, where the genera- 
tors are represented as dominos as follows: 

1 

1 



1 

i+1 



ei = 



D C 



i+1 



(5.2) 

and a product of elements is represented by the concatenation of the corresponding domi- 
nos. Note that we have numbered the left and right "string ends" of the dominos from top 
to bottom, 1 to n. The relation (ii) expresses the locality of the e's, namely that the e's 
commute whenever they involve distant strings. The relations (i) and (iii) read respectively 



O 



= qci = q 



J C 



1 

i+l 



{in) Cj Cj+i ei = 



i+1 



J C 



i+1 



(5.3) 

In the relation (i), the loop has been erased, but affected the weight q. The relation (iii) 
is simply obtained by stretching the {i + 2)-th string. 
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Moreover the algebra is equipped with a trace, defined as follows. Given a domino 
d, we put it on a cylinder by identifying its left and right string ends. We then count 
the number n{d) of connected strings in the resulting picture. The trace is then simply 
Tr((i) = This definition is extended to any element of the algebra by linearity. This 

trace has the important Markov property which allows to compute it by induction: 

Tr(d(ei, 62, ej)ej+i) = ^Tr((i(ei, 62, e^)) (5.4) 

for any algebra element d{ei, Cj) involving only the generators ei, 62, ■■■Cj. 

The crucial remark actually concerns the left ideal l2n{(l) of TL2n{(l) generated by the 
element eie^...e2n-i- The corresponding dominos have their right string ends paired by 
single arches linking the ends 2i — 1 and 2i, i = l,2,...,n. Therefore their 2n left string 
ends are also connected among themselves, thus forming an arch configuration of order 
2n after eliminating the internal loops and stretching the strings (by use of (z) — [in)) 
and giving them a semi-circular shape. We therefore have a bijection between the set 
of arch configurations of order 2n and a basis of the ideal /2n(<?) = ^-^2n('?)eie3...e2n-i5 
made only of "reduced" dominos, namely with all their loops removed and their strings 
stretched. The set of reduced dominos of l2n{,(l) is dnoted by D2n- As a consequence we 
have dim(/2n(i?)) = c„. We denote by the same letter a G or L>2n the arch configuration 
or the corresponding reduced domino. 

5.2. Meander Polynomials 

Given a pair a, h of reduced dominos of l2n{(l), we may form the scalar product 

(a,6) = ^Tr(6*a) (5.5) 

where by 6* we mean the refiected domino wrt one of its vertical edges. The middle of 
the concatenated domino 6*a is nothing but the meander with a as upper-half and 6* as 
lower one (tilted by 90°), while it also contains n small circles formed in the cylinder 
identification, hence 

(a, 6) = Q^^"''') (5.6) 
We deduce the following expressions for the meander and semi-meander polynomials in 



purely algebraic terms ||T2 



m2u{q) = ^ (a, 6) 

aeD2n 



40 



As the trace can be computed by induction, using (p^), we now get an inductive way of 
computing the meander and semi-meander polynomials. This is however far from explicit, 
and does not allow a priori for a good asymptotic study of these polynomials. 

5.3. Meander Determinants 

On the other hand, we may consider other meander-related quantities which are ex- 
plicitly calculable by use of the Temperley-Lieb algebra representation theory, the most 
interesting of them is the "meander determinant" constructed as follows. We first define 
the Gram matrix for the reduced basis of /2n(<?) as the c„ x Cn matrix G2n(<?) with entries 

G2niq)a,b = ia,b) ya,beD2n (5.8) 

The meander determinant is then defined as the Gram determinant 



^2n{q) = det(G2n(g)) 



(5.9) 



This determinant is computed by performing the explicit Gram- Schmidt orthonor- 
malization of the scalar product (., .), and using along the way the representation theory 



of TLn (q) . The result is remarkably simple [0 flS 



A2n(?) = n ^-(5)"^ 



m=l 



2n 
n — m 



2n 

21 J + 

n — m — 1 



2n 

n — m — 2 



(5.10) 



where the Um{q) are the Chebyshev polynomials of the second kind defined by Um{q) 
sin((m -|- 1)6) / sm.[6) with 2 cos 9 = q. 



6. Matrix Model for Meanders 

We have seen in Sect. 2 how to tailor matrix integrals to suit our combinatorial needs. 
The meanders are basically obtained by intersection of two types of curves: the roads 
and the rivers. In the following, we will define a matrix model that generates arbitrary 
configurations of any number of rivers (now viewed as closed curves, i.e. in the case of just 
one river we should close it into a circle), crossed by any number of non-(self)intersecting 
roads lITOf . 
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6.1. The BK^W Model 

Let us now construct a Hermitian matrix integral that generates the meander polyno- 
mials. The computation of such an integral must involve fatgraphs with double-line edges, 
which we will eventually interpret as the river (s) and the road(s). Let us paint in white 
the river edges, and in black the road edges. We therefore have a "black and white" graph 
made of black and white loops which intersect each other through simple intersections. To 
assign a weight say q per black loop (component of road) and p per white loop (component 
of river), the simplest way is to use the replica trick of Sect. 2.3: for positive integer values 
of p and q, introduce q "black" Hermitian matrices Bi, B2, Bq and p "white" Hermitian 
matrices Wi, W2, Wp, all of size N x N, with the only non- vanishing propagators 

white edges : {{Wa)ij{Wb)ki) = j^5a,h5ti^jk ^ 
black edges: {{Ba)ij{Bb)ki) = j;^5a,bhl^jk = \ ^ \ 



(6.1) 



and only simple intersection vertices 



Ay 



Ay 



(6.2) 



for all 1 < X < p, 1 < y < q. 

The case of a unique river will then be recovered by taking the limit p — > 0. This 
suggests to introduce the "black and white" matrix integral 

q p 



-NTiVi{Ba},{Wb}) 



° a=l 6=1 

V{{Ba}dm) = \[Y.Bl + Y.W^ -xY,Y.BaW,BaWk) 



(6.3) 



2 

0=1 6=1 0=1 6=1 



where Zq is a normalization factor ensuring that Zq^p{N:,Qi) = 1. As before, the corre- 
sponding free energy may be formally expanded as a sum over all possible connected black 
and white graphs as 

F<i,p{N\x) = -^LogZq,p{N;x) 

1 (6.4) 

^ |Aut(r)| X (I V 

black &white T ' ^ 
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where, as in Sect. 2.2, g(T),v{T) and Aut(r) denote respectively the genus, number of 
vertices, and symmetry group of F, whereas Lh{r) and Lw{r) denote respectively the 
numbers of black and white loops of resp. black and white edges in F. To get a generating 
function for meander polynomials from ( |6.4|) , we simply have to take ^ oo to only 
retain the planar graphs (genus zero), and then only compute the coefficient of w in the 
resulting expression as a power series of w, which leads to 



N^oo dp 



m2n(<?) (6.5) 



^ 4n 



where we simply have identified planar black and white fatgraphs with meanders (with the 
river closed into a loop) , whose symmetry group is x ^2 for the cyclic symmetry along 
the looplike river, and the symmetry between inside and outside of that loop. The meander 
polynomial is as in (|4.1| ). The large n asymptotic behavior of the meander polynomial 
TO2n(<?) can be directly linked to the critical behavior of the generating function Fq{x) 
i WM - Indeed, (|4.4| ) translates into a singular part 

F,{x),,^^r^{x{q)-xr^'^^ (6.6) 

where 

x{q) = M = a{q) - 1 (6.7) 

So, if we can investigate the critical properties of Fq{x), the meander asymptotics will 
follow. 

The B&W model in the large limit also allows for a simple representation of semi- 
meanders as a correlation function. Indeed, considering the operator 

<Pi = lim ^Tt{Wi) (6.8) 

and computing the correlation function 

(<^i<^i),,P = ^J n W5a</>i</>ie-^^'-^«^^>'{^^» (6.9) 

by use of the Wick theorem, we see that (pi creates a white line at a point, hence the 
fatgraphs contributing to the p — »• limit will simply have a white segment joining the two 
endpoints created by the two ^I's, intersected by arbitrary configurations of road. Given 
such a planar graph, we may always send one of the two endpoints to infinity thus yielding 
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a semi-infinite river, and the fatgraphs with n intersections are just the semi-meanders of 
order n. Hence 

Y^mM^" = {Mi)<ifi (6.10) 

Again, we will get the semi-meander asymptotics from the singular behavior of this corre- 
lation function when x x^- 

Before taking the p — limit, we could have written 



FqA^) = Fq^p{N]x) = V— m2n(?,p) (6.11) 

n=l 

where we have defined a meander polynomial m2n(?,p) for multi-component meanders with 
also multiple rivers (and a weight q per road and p per river). In particular m2n{(liP) ~ 
P"^2n(?) when p — > 0. Similarly, 

(<^i</>i)g,P = mr,{q,p)x'' (6.12) 

n>l 

defines a semi-meander polynomial with also multiple rivers, one of which is semi-infinite. 
In Sect. 6.3 below, we will derive the exact asymptotics of the polynomial Tn2n{Q,P = 1) 
from the Black and White matrix model, while in Sect. 6.4, we will present an argument 
yielding all the critical configuration exponents for large n. 

Note also that if we keep finite, we get an all genus expansion 

aF,,p(A;x) 



dp 



g>0 n=l 



where we have defined the genus g meander polynomials mn\q) that counts the natural 
higher genus generalization of meanders of given genus and number of connected compo- 
nents of road, and with one river. 

6.2. Meander Polynomials and Gaussian Words 

The integral ( |6.3| ) can be considerably simplified by noting that it is just multi- 
Gaussian say in the white matrix vector W = (VFi, Wp). Using the recipees of Sect. 
3.4, we first identify the quadratic form involving each Wb in ( |6.3| ), namely 

N 

W^QWt= Yl iWb)J^Q^JMWb)kl 

hj,k,i=i (6.14) 



= Tt{W^ -xY WbBaWbBa) 



a=l 
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hence 



Qij,kl — ^ik^jl — X '^^{Ba)ik{Ba) 



a=l 



or more compactly 

g 

a=l 

After integration over the VF's in ( |6.3| ), we are left with 
Extracting the p ^ limit is easy, with the result 



(6.15) 



(6.16) 



(6.17) 



dF,^^{N-x) 



dp 



i^(-iTYLogQ) 



(6.18) 



where the bracket stands as in Sect. 2.3 for the multi- Gaussian average wrt the S's. 



Expanding Log Q of ( 6.16 ) as a formal power series of x, we get finally 



9F,,p(iV;x) 



dp 



E 



E 



P—'J m>l ai,a2...,am = i 



(6.19) 



Comparing this with ( |6.13|) , we identify the all-genus meander polynomial generating func- 
tion 

„ 1 



5^iV-2^m(f)(Q)= Yl 



TT{Ba,Ba,...Bajn 



(6.20) 



ai,a2...,a„ = l 

Concentrating on the genus zero case, we may express the meander polynomial (|4.1J ) 
multi- Gaussian integral 



1 I 
m2n{q)= Yl Ji"^ {\-Tr{Ba,Ba,...BaJ\^) 



ai,a2...,a2„ = l 



(6.21) 



where we have noted that the rhs of ( |6.2C1| ) vanishes for large if n is odd, by a simple 
parity argument. Moreover, the rhs of ( |6.2ip must be computed using the planar Wick 
theorem. It is clear however that no pairing can be made between B matrix elements 
pertaining to the two trace terms, as that would violate planarity. In other words, in the 
large N limit, we have {Tr f{Bj)Tr g{Bj)) = {Tr f{Bj)){TT g{Bj)) for any functions / and 
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g of the S's. Hence, using the definition (|3.71|) of Sect. 3.6, we may finaUy express the 
meander polynomial in terms of planar multi- Gaussian averages of words 

1 

m2n{q)= J] |7ai,a2,...,a2„|^ (6.22) 

ai,a2...,a2„ = l 

where we have set for ai, ...Op G {1, 2, q} 




(6.23) 



Fig. 9: A typical pair of Wick pairings used to compute the quantity 
l7i,2,2,i,3,3,i,i,i,i,3,3p- Each pair is represented by a semi-circle with the corre- 
sponding color (l=solid, 2=dashed, 3=dotted). A pairing for 7 is represented, 
as a colored arch configuration of order 12 here. A pairing for 7* is also repre- 
sented, but head-down, to match the bridge colors. The net result is a colored 
meander with 3 components of road. 

Having proved ( |6.22[ ) through some quite lengthy process, it is instructive to interpret 
( |6.22| ) directly. Using the planar Wick theorem to compute the rhs of ( |6.22D , let us represent 
for 7ai,...,a2„ the chain of 2?! matrix elements as equally spaced points along a line, each 
with the corresponding color ai G {1,2,..., q\ (see Fig. 9 for an illustration). A given planar 
Wick pairing is nothing but an arch configuration linking these points by pairs, provided 
they have the same color. As we have to multiply two 7's, we get a pair of such colored 
arch configurations. The fact that the two 7's are conjugate of one another simply means 
that the two arch configurations may be superposed, and that their bridge colors match: 
indeed, 7*^ = 7a2„,...,ai5 so we just have to represent its pairings head down and use 
the same line as for 7. So the net result is the production, for each (pair of) Wick pairing 
pertaining to the two 7's and for each choice of colors ai,a2, ■■■,a2n of the bridges, of a 
colored multi-component meander, namely with colored roads. As the colors are summed 



46 



over, we get a final weight of q per connected component of road, which matches the 
definition of the meander polynomial. 

This interpretation of the result ( |6.22|) allows to immediately generalize it to the 
semi-meander case: 

tT^niq) = ^ ] 7ai,a2,...,a„,a„,a„_i,...,a2,ai (6.24) 

ai,a2,...,a„ = l 

Indeed, the planar Wick pairings for the 7 form again colored arch configurations of order 
2n with the constraint that they have symmetrically identical colors. This latter constraint 
can be represented by a lower rainbow diagram r2n with the corresponding colors. The 
net result is a colored semi-meander, and the result ( |6.24| ) follows exactly like before. 

By analogy with (|6.20| ) , we may also define the genus g semi- meander polynomials as 



J2N-''mi^\q)= E (l^^(Sa,Sa2-Sa„Sa„...SaJ) (6.25) 

g>0 ai,a2...,a„ = l 

Having expressed all our quantities of interest in terms of 7'a, we now have an alter- 
native route for evaluating them. We have to use the recursion relations ( |3.74| ) to first 
compute the 7's and then to substitute them in the above expressions for (semi-) meander 
polynomials. This is still however a tedious work, and it doesn't really allow for any 
asymptotic study at large n. One advantage though of these expressions is that as the 7's 
are positive numbers, they allow for exact lower bounds on the meander polynomials. 

A last remark is in order: in this section, the meander and semi-meander polynomials 
have been computed only for positive integer values of q. Analytic continuation to any 
positive real q is however immediate as these are polynomials. 

6.3. Exact Asymptotics for the case of arbitrary many rivers 



In this section, following we will use the black and white matrix model for mean- 
ders to evaluate the exact asymptotics of slightly different quantities, namely the "multi- 
meander" polynomials m2n(?, 1) counting the meanders with fixed numbers of connected 
components of roads and arbitrary numbers of rivers, closed into non-intersecting loops. 
However this side result already announces the fiavor of what true meander asymptotics 
should look like. 
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We wish to compute the function i ( |6.3|) by first integrating over aU the B matrices, 
rather than the single W one. Using the form (|6.16| ) with B —>■ W, we are then left with 

This is just a Gaussian average over one Hermitian matrix W. Applying the reduction 
method of Sect. 3.1, we arrive at the eigenvalue integral 



This is very similar to the eigenvalue integral ( p.49| ) we obtained in the case of the 0(n) 
model of Sect. 3.4. Actually, upon the change of variables 

1 1-2 

w=^- (6.28) 

we have 

2^(2-9) 

(6 29) 

Let us now evaluate the large A^ behavior of the integral using the saddle-point technique 
of Sect. 3.4. The action to be minimized reads now 

1 ^ 1 

>5(2l,...,2Ar) = — ^^;(2i) - — J2 ^^S\z^-Zj\-^ ^Og{Zi + Zj) 

i=l i<i^j<N l<i,j<N 

v{z) = ^{\^)' + {q-2)Log{l + z) 

(6.30) 

Note that the sum in the last term of 5" may be restricted to i ^ j up to a term of order 
0{1/N) that does not contribute to the large A^ leading asymptotics. Expressing that 
dz^S = 0, we finally get 

This is exactly the saddle-point equation ( |3.51j ) for n = g and the particular choice ( |6.30| ) 
of the potential v{z). Note that v' is a meromorphic function of z with a third order pole 



^2 
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at 2; = — 1. As before, we assume that the hmiting eigenvalue distribution has a support 
made of a single interval [a,6], < a < 6. This requirement turns out to fix entirely the 
meromorphic functions S{z) and P{z) in ( p.62|) . The critical singularity of the genus zero 



free energy / ~ — x) is found to lie at ||25 



X = x{q) 



P 



2 sin(7r|-) 



(6.32) 



where we have set q — 2 cos(7r/), and the corresponding critical exponent 7 reads 

/ 



7 



1-/ 



(6.33) 



As mentioned before, these translate into multi-river meander asymptotics 

R{q, 1)2- 



W2n(?, 1) 



(6.34) 



as 



1 sin^fyr^) 
i?(g,l) = — =2- ^ 2^ 



a(g, 1) 



x{q) 
2-/ 
1-/ 



(6.35) 



In particular, we recover from these values the case of meanders with one river and arbitrary 
many connected components of road, with q = 0, f = 1/2, and i?(l,0) = i?(0, 1) = 4, 
a(l,0) = a(0, 1) = 3. We list a few of these values for various fractions / in Table II 
below. 



Q 


/ 


R{q,l) 


a{q, 1) 


R{q) 





1 
2 


4 


3 


3.50 


1 


1 


9 


5 


4 


3 


2 


2 




V2 


1 

4 


16-872 = 4.68... 


7 

3 


4.13... 


Vs 


1 

6 


36-18^3 = 4.82... 


11 

5 


4.27... 


2 





^ = 4.93... 


2 


4.42... 



Table II: Multi-river meander asymptotics. We have listed a few values of g = 2cos7r/, 
together with the corresponding values of R{q, 1) and a{q, 1), and the numerical values of 
R{q) obtained from the direct enumeration results of Sect. 4.2 for comparison. Note that 
R{q, 1) ~ i?((3') + 1/2 with a good precision. 
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6.4- Exact Meander Asymptotics from Fully-Packed Loop Models coupled to Two- 
dimensional Quantum Gravity 



In a recent work it has been noticed that the B&W matrix model is the natural 
random surface version of some Fully Packed loop model on the square lattice. The latter 
is defined by assigning a color (B or W) to each edge of the square lattice, in such a way 
that two edges of each color meet at each vertex. These edges then form (Fully Packed) 
loops each of which is assigned a weight p or q for W and B loops respectively. The model 
is called the FPL'^{p, q) model [^. When defined on a random surface of genus zero, 
the model assigns colors to the edges of a random fatgraph with only vertices of the form 
BWBW (crossing) or BBWW (avoiding). The B&W model of Sect. 6.1 does not have the 
second kind of vertices. Therefore the original Fully Packed loop model has been further 
restricted. 

The detailed study of the FPL'^{j>,q) model shows two remarkable facts: (i) it is 
critical for all values of < p, < 2 (ii) it is represented in the continuum limit by a 
Conformal Theory with central charge 

cfpl(?, p) = 3 - 6 + (6.36) 

where p = 2 costtc and q = 2 costt/. This was proved by mapping the FPL'^ip, q) model 
onto a three-dimensional height model, where the heights are defined in the center of each 
face, with an Ampere-like rule prescribing the transitions from one face to its neighbors. In 
the continuum limit, the height variable becomes a three-dimensional free field (conformal 
theory with central charge c = 3) , and the corrective weights assigning the factors p and q 
per loop of each color account for the correction of c by electric charges e, /. 

The restriction we impose here on the FPL^{p, q) model on a random surface amounts 
to restricting the height variable to only two dimensions instead of three. The correct 
formula for the flat space Fully Packed Loop theory is therefore 

c(g,p) = 2-6 + 



1-e 1-/ 

p = 2cos7re ^^■^^'> 
q — 2 cos Tvf 

with e, / G [0, 1/2] (i.e. < p,q <2). We therefore state that the B&W model of Sect. 
6.1 is described in the planar (large A^) limit by the gravitational version of a conformal 
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theory with central charge (|6.37|) , namely the same theory defined on fiuctuating surfaces, 
that have to be summed over statistically. Note that the 0{n) model, whose gravitational 
version has been introduced in Sect. 2.3 is also described in the dense phase by a conformal 
theory of central charge c = 1 — 6g'^/{l — g) where n = 2cosng. So we are now dealing 
with a sort of double 0{n) model coupled to gravity. 

The coupling to gravity of a conformal theory with central charge c < 1 has been 
extensively studied within the context of non-critical string theory. The gravitational 
theory has a new parameter x, called the cosmological constant, coupled to the area of 
the surfaces we have to sum over. More precisely, the free energy for a conformal theory 
coupled to gravity in genus zero reads 

F-LogZ=^x^ J2 ZcFTir) (6.38) 

^ Q connected surfaces F 

of area A 

where Zcft{^) denotes the partition function of the conformal theory on the genus zero 
connected surface F. Comparing Z with the B&W model partition function ( |6.3|) , we see 
that X plays the role of cosmological constant, as n = A are the areas of the tessellations 
dual to the fatgraphs of the model. When the conformal theory has central charge c, the 
free energy ( |6.3(j| ) has been shown to have a singularity of the form ( |3.63| ) ( P^BBf ): 



Fr^ixc-xf--^ 7 = ^(c-1- V(l-c)(25-c)) (6.39) 

when X approaches some critical value Xc- This is easily translated into the large area 
asymptotics of the partition function of the model on surfaces of fixed area 

F^-^ (6-40) 

Moreover, the operators of the conformal theory get "dressed" by gravity, and their corre- 
lation functions have singularities of the form 



{(Pm,...(pm,) ~ {Xc - x)E^-.-^+2-'= (6.41) 



where the "dressed dimensions" are related to the conformal dimensions hm of their 
undressed versions in the conformal theory through 

A. = (6.42) 
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We have now all the necessary material to compute the configuration exponents of all 
the meandric numbers of interest. Applying the result ( |6.4U| ) to the central charge ( |6.37| ), 
we find the configuration exponent of the multi-river meander polynomial ( |6.11|) 

6.43 



a{q,p) = 2 + ^v^l - c(q,p)(v/25 - c{q,p) + - c{q,p)) 



as well as that of the multi-river semi-meander polynomial (|6.12 

_ Rjq.pr 

a{q,p) = a{q,p) - 1 + 2Ai 

where Ai is the dressed dimension of the operator creating a white endpoint. In the 
conformal theory, this operator is known to have the dimension 

= ^ - 4(1^ 

where p = 2 cos ne. We simply have to apply (|6.42|) to (|6.45|) and substitute the value of 
Ai back into j ^M )- 

For p = 1 (e = 1/3) and q arbitrary, we find 

2 - / 

^ (6.46) 



a{q, 1) 



1-/ 

for all q = 2 cos it f. The first line of ( |6.46|) agrees with the saddle point result ( |6.35|) . 
The second line is readily obtained by noticing that Ai = when e = 1/3, and therefore 
a{q, 1) = «((?, 1) - 1. 

For p = q' = 0(e = / = l/2),we get the exact values of the meander and semi-meander 
configuration exponents: 



a(0,0) =2 + -^V5(V5 + v^) 



Y (6.47) 
a{0,0) = 1 + —Vn{V5 + V29) 

Note that the arguments of this section do not give any prediction for non-universal 
quantities such as R{q, p) (which is expected to depend on q and p explicitly, not just on 
c{q,p)). 
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Finally, for p = and q arbitrary, we find: 

a{q, 0) = 2 + ^V1-c(q)(V25-c(q) + v^l - c{q)^ 
a{q, 0) = 1 + ^V3-4c(q)(^V25-c(q) + v^l - c(q)^ 



with c(q') given by ( |3.69D . Note that the second line of ( |6.48D breaks down when q = qc 
corresponding to c{qc) = 3/4, namely 



/ ,/97 — 1 \ 

gc = 2cos(7r ^ j (6.49) 

We identify this as the critical value of q beyond which the winding becomes relevant in 
semi-meanders, namely when circles (roads intersecting theriver only once) dominate the 
semi-meander configurations. 



53 



PART C: Fluid Membrane Folding 



We now turn to our second application of matrix integrals, having to do with the 
generation of foldable two-dimensional triangulations. These triangulations are a simple 
model for so-called tethered or fluid membranes, objects of physical and biological interest. 

Although we now deal with the folding of two-dimensional objects (as opposed to the 
one-dimensional polymers of part B), we will rephrase the problem as that of enumerating 
vertex-tricolored triangulations. A suitable matrix model will be presented and solved 
using various techniques. 



7. Folding Triangulations 

As mentioned above, one usually distinguishes between two distinct models for mem- 
branes: 

Tethered Membranes: are represented by reticulated networks, typically some domain 
of a regular lattice with physical vertices and rigid bonds, allowed only to change their 
spatial configuration through folding. 

Fluid Membranes: are represented also by now irregular reticulated networks in which 
the valency of vertices is no longer fixed, this disorder accounting for the fluidity. But the 
bonds are still rigid, and the membrane is still allowed to change its spatial configuration 
through folding. 



In the following, we first briefiy describe some results on the triangular lattice 
folding (tethered membrane folding), before turning to the problem of folding triangula- 
tions (fiuid membranes). Here we only consider the so-called "phantom folding" of 
membranes, i.e. we allow the network to interpenetrate itself, and are only interested in 
the statistics of the final folded states, independently of the actual feasability of the folding 
process. Moreover, the folding of the triangular lattice will be two-dimensional, in that 
the folded configurations will be subsets of the original lattice. 
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(a) 



(b) 



Fig. 10: A choice (a) for tiie tangent vectors of tlie triangular lattice, to- 
gether with the corresponding coloring of the edges by 1, 2, 3. This is the flat 
configuration of the membrane. A folding configuration (b) with the corre- 
sponding folded bonds (thick black lines) and edge coloring. The three colors 
correspond to the three unit vectors with vanishing sum represented above. 

7.1. Folding the Triangular Lattice 

We consider the triangular lattice, or rather a rhombus-shaped portion of it with size 
N X P. Let us view this as a reticulated network, with solid edges of unit length. Each 
such edge carries a unit tangent vector t, with the constraint that 



With the choice of tangent vectors indicated in Fig. 10 (a), a folding of the triangular lattice 



unit vector el, we see that the images of the tangent vectors of S may only take three 
values el, 62, 63, where el, 6'2, es are three unit vectors with vanishing sum, hence forming 
angles of 27r/3. Let us associate colors numbered 1, 2, 3 to these three possible images. 

A folding map p of the triangular lattice is therefore a coloring of its edges, with the 
three colors 1, 2, 3, such that the three colors of edges around each face are all distinct. An 
example of such a coloring is given in Fig. 10 (b) together with the corresponding folding 
configuration. The dual of this coloring model is the problem of tri-coloring the edges of the 
hexagonal (honeycomb) lattice in such a way that the three edges adjacent to each vertex 
are painted with distinct colors 1,2,3. It has been solved by Baxter |]3^, by use of the 
Bethe Ansatz. Baxter's results yield in particular the exact value for the thermodynamic 
entropy of folding per face of the triangular lattice s = limAr,p^oo j^^og Zn^p, where 




(7.1) 



around faces 



is a continuous map p : S ^ IR^, preserving the length of the tangent vectors and satisfying 
the condition ( [7.1|) around each triangular face of S. Let ti,t2,i3 denote the unit tangent 
vectors to a given face of S. Their images p{ti) are three unit vectors with vanishing sum, 
according to the face rule ( |7.1| ). Fixing the image of one tangent vector of 5" to be a given 
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Zn^p denotes the total number of distinct folded configurations of our portion of lattice. 
It reads 

s = Log(|^r(l/3)^/2) (7.2) 

This was originally proved by explicitly diagonalizing a large (transfer) matrix, indexed 
by the coloring configurations of rows of edges in the honeycomb lattice, and describing 
the "row-to-row transfer" , i.e. the allowed coloring configurations for two such neighboring 



rows. The thermodynamic entropy (7^) is then the logarithm of the largest (Perron- 
Frobenius) eigenvalue of this matrix. The diagonalization is performed using a particular 
ansatz for the form of the eigenvectors, the Bethe Ansatz. The proof of (|7.2|) being highly 



technical, we will not reproduce it here, but refer the interested reader to the original 



papers |3y]. Let us simply mention that this model is part of the class of Two-dimensional 
Integrable Lattice Models, for which a Bethe Ansatz solution exists. 

A score of other lattice folding problems have been studied [ETl. Actually, one can 



classify all the (compactly) foldable lattices ||3^, and define higher-dimensional generaliza- 



tions of their folding problems. Remarkably, all of them lead to some equivalent coloring 
problems, also rephrased into Fully Packed loop enumeration problems. For instance, in 
the case of folding the triangular lattice, we have seen that the model is equivalent to that 
of tricoloring the edges of the lattice. Concentrating say on the colors 1 and 2, we see 
that edges of alternating colors 1212... form loops on the dual (honeycomb) lattice, and 
moreover each vertex of the dual is visited by exactly one of these loops: the loops 1212... 
are therefore fully packed, as well as the 2323... and 1313... Remarkably, the same type of 
Fully Packed loop models have emerged in our study of meanders (see Sect. 6.4 above). 

7.2. Foldable Triangulations 

Fluid membranes are modelled by irregular networks of vertices linked by edges, in 
which the valencies of the vertices are arbitrary, as well as the genus of the underlying 
surface, which might have an arbitrary topology. As advocated in Sect. 2, matrix models 
provide us with a means of generating such graphs. 

Let us restrict ourselves to networks with only triangular faces, namely triangulations. 
The initial question one can ask is: are all triangulations foldable? It is indeed desirable 
to consider only triangulations with a large number of folded configurations, otherwise the 
effect of folding might be wiped out in the limit of large size. We therefore demand that 
our triangulations be compactly foldable, namely that one can fold them completely onto 
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just one of their (equilateral triangular) faces. With this constraint, it is clear that not 
all triangulations turn out to be "foldable". Indeed, let us paint by three distinct colors 
1,2,3 the three vertices of the image triangle, and paint accordingly the vertices of the 
preimages under the folding map. This results in the tri-coloring of the vertices of the initial 
triangulation, in such a way that the three colors around each triangular face are distinct. 
So only the vertex-tri-colorable triangulations will be compactly foldable. Another way of 
viewing this restriction is to recall that we first need to attach tangent vectors to the edges 
of the triangulation, in such a way that ( [7.1|) is satisfied. It is straightforward to see that 
this is possible only if the vertices of the triangulation are all even, as around such a vertex, 
we must have an alternance of tangent vectors pointing to and from it. This condition 
turns out to be sufficient in genus zero to grant the tri-colorability of the triangulation. 



The situation in higher genus is unclear [33 



Let us now introduce a generating function Z{xi, x^-, X3; t; A^) for possibly disconnected 
vertex-tricolored triangulations of arbitrary genus, such that 

F(xi, X2, xa; t; A^) = Log Z(xi, X2, X3; t; A^) 



vertex — tricolorcd 
connected triangulations T 



^1 ^2 ^3 - |Aut(T)| 



^ ^.,(T)^..(T)^.3(T)t^iV'-''^(^) (7.3) 



where ni{T) denote the total numbers of vertices of color z, A{T) the total number of 
faces, h{T) the genus and |Aut(T)| the order of the symmetry group of the tricolored 
triangulation T. 

The construction of a matrix model to represent Z'(xi, X2, X3; t; A^) is based on the 
following simple remark: in a given tri-colored triangulation T if we remove say all the 
vertices of color 3 and all edges connected to them, we end up with a bi-colored graph, with 
unconstrained vertex valencies. Such bi colored graphs are easily built out of the Feynman 
graphs of a two Hermitian matrix model, say Mi and M2, the index 1 and 2 standing for 
the color, with colored vertices 



V*, *" <<>>^ 



57 



(7.4) 



connected through propagators imposing the alternance of colors, namely 
{{Ma)^j{Mh)ki) = {l-Sab)SjkSu^ = :=^-::::l 

(7.5) 

where a,b = 1,2. Let us introduce the corresponding matrix integral, but keep fixed 
while the matrices are taken of size n x n, n possibly different from N. This gives the 
partition function 

y(Mi,M2;xi,X2,t) = xiLog(l - Ml) + X2Log(l - M2) + ^MiM2 

where the normalization factor (f>n{t,N) ensures that Z^(0,0;t; A) = 1. Clearly we must 
take the integral ( [7.6|) only at the level of formal series, by expanding it in power series 
of xi and X2, and computing the Gaussian integrals with measure dMidM2e~^^^^^ by 
the Feynman procedure. In particular, this gives the following rule for two-dimensional 
integrals 

= — f dxdye'^'^y x'^y'^ 

/ ^ X a (7.7) 
= 5«^r(a + l) ^- 

The integral could be made rigorous and convergent by considering instead a normal matrix 
Ml — > M, M2 M*, but we will content ourselves with formal power series anyway. 
Let us now compute the Feynman graph expansion of the free energy 

Fn{xi,X2]t]N) = LogZn{xi,X2;t;N) 

y -.(r)ii^(r)^(i^(r)-i.(r))^F(r) (7.8) 

, ^ Aut F ^ 2 

bicolored cormected ' \ / 1 

graphs r 

where we have denoted by ni(F) the number of vertices of color z, V{r) = ni(F) + n.2(F), 
E(r) the number of edges, -F(F) the number of faces (the boundary of which are loops of 
running indices, accounting for a factor n each), and |Aut(F)| the order of the symmetry 
group of the bicolored fat graph F. Adding a central vertex of color 3 in the middle of 
each face of F, and connecting it to all the vertices around the face with edges will result 
in a vertex-tricolored triangulation T. The number of such added vertices is nothing but 
ns{T) = F{T). Introducing 

X3 = ^ (7.9) 
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we may rewrite 

Fr,{xi,X2;t;N) = F{xi,X2:X3;t;N) (7.10) 

by use of the Euler relation 2 - 2/1(7) = 2 - 2h{T) = V{r) - E{r) + F(r) and the fact that 
A{T) = 2E{r), as each edge of F gives rise to two triangles of T, one in each of the two 
faces adjacent to the edge. It is also a simple exercise to show that |Aut(T)| = |Aut(r)|. 
Considering ( |7.10| ) as a formal power series of t with polynomial coefficients in xi, 0:2, X3, we 
see that the knowledge of -F„(xi, X2; t; N) for integer values of n determines completely the 
polynomial dependence on X3 = n/N at each order in t. Hence computing F„(xi, X2; t; N) 
through the integral formulation (|7.6| ) will yield the generating function for compactly 
foldable triangulations. 



8. Exact Solution 

We now present the exact computation of the generating function F{xi,X2:Xs;t; N) 
for compactly foldable triangulations. Although it can be obtained through orthogonal 
polynomial techniques generalizing that of Sect. 3.2, we choose to present an alternative 
and powerful approach, using the Discrete Hirota equation, in the form of a recursion 
relation for the quantities Fn{xi,X2;t] N), in terms of n,Nxi and Nx2- We will then 
write some direct formal series expansion for the solution, and use saddle-point techniques 
to extract the large N behavior of the free energy. This will give a nice formula for the 
generating function of vertex-tricolored triangulations of genus zero. 

8.1. Discrete Hirota Equation 

The step zero in computing the integral ( [7.6|) is like in Sect. 3.1 the reduction to eigen- 
value integrals. Diagonalizing both matrices as Mi = UirriiUj , mi diagonal, Ui unitary, the 
change of variables (Mi,M2) (C/i, mi; t/25 ^^2) has the Jacobian J = A(mi)^A(m2)^. 
We encounter a problem as the potential V{Mi, M2:, xi, X2,t) is not invariant under uni- 
tary conjugation, namely the term Tr(MiM2) = Tr(0mi0^m2), with the unitary matrix 
Q = U2U1. So the integral over Ui is no longer trivial, and yields a factor 

/c/ne-^T^(""^^"'"^^» oc det(e-f-^^'-^.Oi<.,,<. 
J A(mi)A(m2) ^^"'^ 

up to some normalization factor depending on n and only. Eqn. ( pT|) is the celebrated 
Itzykson-Zuber formula of integration over the unitary group [^5|, itself a particular case 
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of the Duistermaat-Heckmann localization formula [36], as the determinant in ( |8.1|) may 



the 



be viewed as a sum over the critical (saddle-) points of the integrand, at 0, 
permutations of the eigenvalues. 

The integral ( |7.6|) then reduces to 



Z^(xi,X2;t;iV) = -^—r / (imidm2A(mi) A(m2)e-^^'-(^("^i'"^^'"i'"^'*)) (8.2 

'^n{t, N) J 



where the normalization factor i/'n(t, n) ensures that Zn{0, 0; t; N) = 1. It is easily derived 
by expanding the two determinants as sums over permutations, with the result 



n „ 

sgn(ar)]^ dmi^idm2,irnl^-^~^ml''^l^~^e~^'^^-' 



7712, 



H-i)l 



n 



n 



L {N/ty 



by use of (fTTI). 
Writing 



A(mi) A(m2) 



M.N) 

n 

X J]e-^"^i-"^^-(l - mi,,)-"(l - m2,^)-'dml,dm2,^ 



1=1 



with a = Nxi, b = Nx2, and using the basic definition of determinants 



Y[{l-mk,^)-'"'A{l-mk) = det 



i=l 



(1 - mk,i) 



j-afe-l 



- l<i,i<n 



J2 sgn(a)n(l-^M)'^^^^""'="' 

O-ESn 



2 = 1 



for k = 1,2, and the shorthand notation 



(8.3) 



•4) 



(8.5) 



Zn{a,b) = Zn{xi,X2;t; N) a = Nxi,b = Nx2 



(8.6) 
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we finally get 



1 /■ 



(1 - mi,i)'^(^)-"-^(l - m2,i)^(')"^"^e-^'^i-'""2- 



n! V- , , (8-7) 
2^ sgn(z/) 



V'n(t,iV) ^^^^ 



n „ 



v{i)-b-l^-^xy 



where we have set v = Ta~^ , with the same signature as err, and explicitly factored out 
the sum over a. Moreover, the dummy integration variables have been rebaptized x and 
y, and the integral can be computed by expanding the integrand as a power series of x, y 
and then using term by term the prescription (|7.7| ) . The partition function takes therefore 
the form 

Zn{a,b) = ■ D^{a,b) (8.8) 

Vn{t, N) 



where Dn{a, b) is the n x n determinant 



Dn{a,b) = det 



j dxdy{l - xy-^-^l - y) 



j-b-l^-Nuxy 



(8.9) 



and Vn(t, N)/n\ = 15.(0, 0) = ni<.<n(^ - ^V-fiN/ty. 

The Hirota equation is simply the rephrasing in terms of Dn{a, b) of the following 
identity satisfied by any (n + 1) x (n + 1) determinant D and the minors Dij obtained by 
erasing the i-th row and j-th column, as well as the minors -Di^,i2;ji,j2 obtained by erasing 
the rows ii, 12 and columns ji, in D: 

-D -Dl,n+l;l,n+l = -Dl^l — -Dl,ri+1 -Dn+1,1 (8.10) 

When expressed in terms of D = Dn+i{a + 1,6+ 1), this implies the quadratic equation 

Dr,+i{a + l,b + l)D^-i{a,b) 

(8.11) 

= L>„(a + 1, 6 + l)I5n(a, b) - D^{a, b + l)Dr,{a + 1, b) 

Finally, using ipn+iit, N)ipn-i{t: N)/iJ;,^{t, N)'^ = n/{Nu) = nt/N from (^), we get the 
Hirota Bilinear equation for Zn{a,b) by substituting (|8.8|) into ( |8.11J ): 

Z^+i(a + l, 6+1) Z„-i(a, b) = Z^(a + 1, 6+1) Z^(a, b)-Z^{a, 6+1) Z^(a + 1, 6) (8.12) 
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This recursion relation determines the partition function Zn{a,b) entirely, once we apply 
the following initial conditions. Let 

Fr,{a,b) = Log Zr,{a,b) = J] (t/iV)-cc;^(a, 6, n) (8.13) 

m>l 

be the generating function for connected tricolored triangulations, with a = Nxi, b = Nx2 
and n = Nxs- Then, introducing the shorthand notation 6xf{x) = f{x + l) — f{x) for finite 
differences, the Hirota equation ( |8.12D turns into a finite difference equation for Fn{a,b): 

SaSMa,b) = -Log(l-n^e^"^"("+^'^+^)-^-^"-i("'^)) (8.14) 

This turns into a nonlinear recursion relation for the coefficients ujjn{o,^b^n). But 
thanks to their interpretation in terms of tricolored triangulation counting, namely 
Y.uj,n{NxuNx2,Nx:i){t/N)'^ = Y.N'^~'^^x1^xl^x';''i^{lr\co\.lriajvg.), all of these are 
polynomials of a, 6, n, and all have at least a term abn in factor, as there is always at 
least one vertex of each color 1,2,3 in such a triangulation. This allows for performing 
the discrete integration step involved in the recursion, and yields the complete solution for 
-F„(a, b) as a formal power series of t. The first few terms read: 

ui = nab 

nab . , , 

UJ2 = -^(n + a + 6) 

= ^^^n" + ^{a + b)n + a^ + ^ab + b'^ + 1) (8.15) 
3 

T}j(xt) 

= — (n^ + 6(a + 6)n^ + (6a^ + 17a6 + 66^ + 5)n 
+ (a + 6) (a^ + hab + 6^ + 5)) 

Note the symmetry in a, 6, n, now manifest. In the limit of large N , the genus zero free 
energy 

fo{xi,X2,X3]t) = lim -^FAr^3(A^xi,iVx2) (8.16) 

N^oo iV 

satisfies the differential equation 

9.,a,j0 = -Log (1-^36^^3(9^1 +9^2 +9.3 )/o) (8.17) 



as a consequence of (|8.14| ) 
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8.2. Direct Expansion and Large N Asymptotics 

Let us explicitly compute the partition function Zn{a,h) (|8.8| ) by expressing the de- 
terminant (|8.9| ) as a formal series expansion in t. In the integrand of ( |8.9| ), we expand all 
terms of the form 

A;>0 ^ ' 

for (w, c) = (x, a + 1 — i) and (y, 6 + 1 — j) respectively, and then compute the defining 
integral by use of the prescription ( [7.7| ). After a little algebra, we arrive at 



Dn{a,b) = det 



^ 1 r{k + a-i + l)r{k + b- j + 1) . t .k+i 
'-k>o ^ ' \ ' JJ 



V- ^ (t/N)^'^+^T(k, + a-i + l) ^ , 



^.19) 



fci ,...,A;„>0 '1=1 



Factoring r(ki + b — n + 1) out of each line of the remaining determinant, we are left with 
the computation of the determinant 

det[{k, + b-j){h + b-j-l)...{k, + b-n + l)] = det[qn-j{h)] (8.20) 

where the polynomials Qm are monic of degree m, and as such satisfy det[qn-j{ki)] = 
det[kj~'^] = A{k) from ( |3.4|) , the Vandermonde determinant of the /c's. Repeating the 
above trick with the columns of the determinant, we finally get the symmetric expression 

Z(ab)- V ±(h)' TT (V^)'^^'"^ + g - n + 1) m + b-n + l) 

A;i,...,fc„>0 1=1 

By construction, this formal series of t is a solution to the Hirota equation ( |8.12| ). Note 
the remarkable similarity between the expansion ( |8.21|) and the reduction to eigenvalues 
of the one-matrix integral ( p.3|) , except that the integration over eigenvalues is replaced 
by a sum over non-negative integers. We indeed have the same two types of terms: (i) 
the squared Vandermonde determinant, i.e. the repulsion term, and (ii) the product of 
ratios of gamma functions, that can be exponentiated into a potential term. The expansion 
( fj.21| ) allows us for a computation of the large limit of the free energy, by applying the 



saddle-point techniques. Indeed, writing ki = Nai, we may estimate the expansion by 
a real integral over the ct's when N is large, of the form J dai...danex.\){—N'^S{{ai})), 
where the index i itself ranging from 1 to n = Nx2, may be written as z = Ns, s G [0, 1]. 
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Expressing that this integral is dominated by the minimum of S, it is a straightforward, 
though tedious calculation to get the large asymptotics of Z. We refer the reader to 
ref. |2^] for details, and simply give the beautifully simple result here. The genus zero free 



energy of the model ( ^.161 ) satisfies the following 

(8.22) 





-F2 


-F3) = txi 


F2{1 


-Fs 


-Fi) = tX2 






-F2) = tX3 



where Fi = Fi{xi, X2, X3; t) = tXi + 0(t^) are formal series of t with polynomial coefficients 
of the x's. Actually it is easy to see that Fi is the generating function for tricolored rooted 
trees, whose root has color i. This result still awaits a good combinatorial interpretation. 
One way of proving the result ( |8.22| ) a posteriori is to check that it satisfies the differential 
equation ( |8.17| ). As an independent check, in the particular case xi = X2 = X3 = z we find 
Fi = F2 = F3 = F and F(l — 2F) = tz, in agreement with a former result of Tutte . 



Analyzing the critical properties of /o as a function of t, we find a singular behavior 
of the form fo{t) ~ {tc — t)^/^ for generic values of the x's {xi > and Xi 7^ Xj for i 7^ j). 
This leads to a number of tricolored triangulations behaving as T(xi, 0:2, xa)^/^^/^ in 
terms of their fixed area A. We recover therefore the value 7 = —1/2 ( |3.68| ) of the pure 
gravity, meaning that the tricoloring constraint does not affect the configuration exponent 
of triangulations, but does affect the leading behavior through the function T(xi, X2, X3) 7^ 
12. 



9. Conclusion 

The combinatorial applications of matrix integrals are various and many, and it would 
be impossible to even enumerate them all here. In these notes we have chosen to concen- 
trate on folding problems with simple physical interpretations, but this is mainly a matter 
of taste. 

In the study of meanders, we have shown that even an almost one-dimensional problem 
had to be first viewed as a random graph (or surface) problem, and then matrix and 
quantum gravity techniques have allowed for obtaining many interesting results. 
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In the case of foldable triangulation enumeration, the matrix model has allowed for 
the derivation of a very simple formula for the genus zero counting function, that still has 
to be interpreted combinatorially. 

The main lesson we would like to draw is that the matrix models give a maybe less 
intuitive but definitely different angle on graph-related combinatorial problems. 
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